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ABSTRACT 


Loran-C  time  difference  (TD)  readings  were  taken 
consecutively  at  several  survey  stations  near  Monterey,  CA, 
over  a  period  of  three  days.  One  station,  Range-7,  was 
designated  the  "known"  point  of  a  differential  Loran-C 
system.  Readings  from  the  known  point  were  used  to  correct 
readings  from  a  second  survey  station.  Luces  Point.  A  method 
of  improving  the  precision  of  Loran-C  TD  readings  based  on 
the  redundancy  and  relative  accuracy  of  three  LOP's  was 
developed  and  applied  to  the  data.  Since  only  one  receiver 
was  available,  a  linear  regression  of  TD  vs  time  was 
calculated  and  used  in  the  differential  correction. 

Based  on  496  sets  of  data  taken  at  6-second  intervals  at 
Range-7  in  two  groups,  before  and  after  about  250  readings  at 
Luces  Point,  the  absolute  accuracy  of  Luces  Point  data  was 
improved  from  about  385  m  to  about  48  m  '  pared  to  the  known 
position  of  the  point.  Precision  was  proved  from  about 
14.9  m  to  about  12.6  m  circular  error  (CEP)  using  the  three- 
station  correction.  Further  improvement  would  probably  have 

resulted  if  two  receivers  were  available  for  real  time 
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I.  nmooucncN 


A.  Purpose 

Hie  purpose  of  this  study  was  to  investigate  the  potential 
acxaoracy  of  the  Loran-C  navigation  system  using  a  differential 
cxirrection  technicjue  and  a  correction  for  the  use  of  three  secondary 
stations  rather  than  the  usual  two. 

Sources  of  error  were  first  examined  cind  their  magnitudes 
estimated.  Methods  of  expressing  the  accuracy  of  a  determination  of 
geographic  position  were  also  considered. 

It  was  determined  that  position  accuracy  is  composed  of  two 
factors:  the  absolute  accuracy  of  the  position  (which  can  be 
considered  the  bias,  or  systematic  error,  in  the  system)  and  the 
precision  of  the  individual  readings. 

Over  longer  periods  of  data  collection  it  was  considered  pxissible 
to  subtract  out  the  system  bias  by  means  of  a  quasi-differential 
correction  based  on  intermittent  time  difference  readings  at  a  fixed 
point,  designated  the  "kncwn"  point  of  the  differential  system,  and 
to  inprove  the  precision  of  TD  readings  by  taking  advantage  of  the 
redundancy  of  the  reception  of  three  master-secondary  pairs. 

For  this  purpose,  a  theory  of  corrections  was  developed  and  cin 
ejperiment  was  designed  to  determine  the  effectiveness  of  the 
corrections.  The  experimental  data  was  processed  by  ocnputer,  and  the 
relative  acxruracies  were  calculated.  Computer  programs  for  processing 
and  displaying  data  were  written  for  this  purpose. 
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B.  Method 


To  collect  data  for  processing,  several  known  survey  points  with 
published  positions  were  chosen  as  locations  for  collecting  Loran-C 
data.  Using  known  points  provided  a  means  of  comparing  results.  One 
point  in  a  relatively  central  location  was  chosen  as  the  "known"  point, 
where  a  larger  amount  of  data  was  collected.  Data  from  the  other  points 
was  corrected  based  on  the  data  from  the  "known"  point. 

During  periods  of  about  ei^t  hours,  the  receiver  was  moved  from 
one  survey  point  to  ainother.  Ihe  emtenna  was  set  up  over  the  points 
and  time  difference  data  was  collected  for  periods  of  about  20  to  30 
minutes  at  each  point.  The  time  required  to  move  the  receiver  frcm  one 
point  to  another  Vcuried  from  about  30  to  60  minutes. 

Data  was  collected  by  a  Racal  Megapulse  monitor  receiver  with 
internal  clock  accuracy  of  about  0.02  jusec.  A  Silent  700  data  terminal 
was  used  to  collect  the  data  on  cassette  tapes.  Data  was  recorded  at 
intervals  of  5  seconds.  Later,  the  data  was  transferred  to  5V'  disks 
for  use  on  a  oorputer. 

Since  there  was  no  way  to  collect  time  difference  data  at  two 
points  simultaneously,  a  method  of  estimating  the  time  differences  at 
the  knavn  point  was  developed.  This  involved  using  the  available  data 
to  establish  linear  regression  lines  of  time  difference  versus  time  at 
the  known  point. 

Using  the  known  time  of  data  collection  at  the  other  points,  an 
estimate  of  the  time  differences  at  the  known  point  could  be  used  to 
calculate  a  correction  for  data  at  the  other  points  at  the  times  of 
data  colleotion. 
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A  number  of  cxsrputer  programs  were  written  to  calaalate  the 
required  information  for  differential  corrections.  A  Correction  For 
Three  Secondary  Stations  was  developed  mathematically  and  used  to 
process  the  data  in  cin  effort  to  reduce  the  random  error  of  the  time 
difference  readings. 

Based  on  reduction  of  standard  deviations  of  the  data  and 
reduced  circular  error,  the  three  station  correction  was 
effective  to  some  degree,  although  absolute  error  turned 
out  to  be  the  limiting  factor  in  the  experiment.  Reduction 
of  absolute  error  by  means  of  the  differential  correction 
was  significant. 
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n.  IHE  ICRAN-C  NAVIGAnCN  SYSEEM 


A.  Basic  Descripticn 

loran-C  is  an  electronic  positicaiing  system  based  on  the  reo^jtion 
of  pulsed  If*  (90-110  }dfe)  signeils  from  a  master  transmitting  station 
cind  one  or  more  secondary  stations  vAiich  ocrprise  a  Loran-C  "cdiain''.  A 
brief  description  of  the  system  will  be  provided  here  as  a  basis  for 
later  discussions  of  the  accuracy  considerations  and  differential 
methods  used  in  the  experiment. 

In  the  hyperbolic  mode,  geographic  position  (GP)  is  determined  from 
time  difference  (TD)  readings  provided  by  a  Lorein-C  receiver.  Ihe 
difference  in  rec^jtion  time  of  pulses  received  fran  the  master  station 
and  pulses  received  frcsn  each  secondary  station  give  the  receiver's 
position  in  terms  of  a  hyperbolic  lattice  si’stQm.  Geogra^iic  position 
can  be  calculated  fron  the  TD  readings  using  loran-C  lattice  tables 
published  by  the  Defense  Mapping  Agency  or  corputed  directly  using 
various  ccnputer  techniques.  Nautical  charts  overprinted  with  loran-C 
TD  lattices  eire  cilso  available  for  plotting  geographic  position 
directly  without  converting  first  to  latitude  eind  longitude. 

Ihe  hyperbolic  lattice  system  is  a  series  of  hyperbolas.  A  pulse 
(actucilly  a  group  of  pulses)  is  first  transmitted  by  the  master 
station.  Ihe  wave  travels  to  each  secondary  station  of  the  loran-C 
chain,  and  each  secondary  transmits  its  cwn  pulse  after  a  time  delay 
(called  the  coding  delay)  vhich  is  unique  to  each  secondary. 
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Ihe  master  pulse  is  received  first  by  the  receiver.  A  short  time 
later  the  pulses  from  the  seoondeuries  are  received.  The  time  dif¬ 
ference  between  reception  of  the  master  and  each  secondary  is  measured 
in  the  receiver.  Each  TD  identifies  one  line  of  position  (IDP)  in  the 
set  of  LDP's  for  a  iraster-secondary  pair.  Figure  II-l  shews  a  cypical 
loran-C  master-secondairy  pair  with  a  set  of  hyperbolic  LDP's.  Each  LDP 
in  this  set  represents  a  line  of  equal  time  difference  between 
reception  of  the  master  pulse  cind  reaction  of  the  pulse  frem  the 
secondary.  Figure  LI-2  shews  an  exaitple  of  a  hyperbolic  lattice  system 
formed  by  two  master-secondary  pairs 

As  with  latitude  and  longitude,  two  LDP's  eure  required  to  define  a 
geographic  position.  Unlike  latitude  cind  longitude,  two  TD  readings 
from  different  secondaries  do  not  uniquely  define  a  GP,  since  the  TD's 
may  intersect,  at  two  points  within  the  coverage  curea.  Ihis  is  normally 
not  a  prctolem,  since  even  an  approximate  kncwledge  of  the  geographic 
position  of  the  receiver  will  cLLlcw  determination  of  the  correct 
intersection  point. 

When  more  than  two  LDP's  are  available  to  define  a  GP  (when  signals 
from  more  than  two  secondary  stations  can  be  received) ,  the 
intersection  of  the  three  LDP's  define  three  geograpjhic  points  (s'*"' 
Figure  IV-2) .  Ihe  accuracy  of  the  two  best  pairs  of  LDP's  is 
sufficient  for  most  navigation  pjurposes,  but  it  is  possible  to  combine 
three  or  more  LDP's  to  estimate  the  position  of  one  point.  A  method 
using  corrections  frem  three  LDP's  is  used  in  the  experiment  and  will 
be  discussed  in  Section  IV. B. 
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B.  Time  Difference  Bquaticxis 

The  difference  between  arrival  time  of  the  signal  frcm  the  master 
and  arrival  time  of  the  signal  from  a  seocxxJary  can  be  described  by 
the  following  equaticxis  (Pef .  1) : 

The  basic  time  'difference  equaticsi  for  Loran-C  is: 

TO  =  ts-t^  [II-l] 

where: 

TD  =  time  difference  at  a  given  point  P. 
tg  =  time  of  reception  of  secondary  pulse, 
t^  =  time  of  rec^jtion  of  master  pulse. 

A  mere  useful  form  of  this  equation  takes  into  account  the 
corponents  which  make  vp  the  travel  time  of  the  signals  over  each 
section  of  the  paths  from  the  stations  to  point  P: 

TD  =  (Tg  +  OTg  +  ED)  -  (Tin  ^^m)  +  [II-2] 

where: 

TD  =  Time  Difference  at  a  given  point  P. 

Tg  =  Distance  in  time  units  between  the  secondary  station 
and  point  P. 

Tm  =  Distance  in  time  units  between  the  master  station 
and  point  P. 

CTg  =  Secondary  Phase  Correction  for  all  seawater  path 
between  secondary  station  and  point  P. 

DTjn  =  Secondary  Phase  Correct iwi  for  all  seawater  path 
between  master  station  and  point  P. 

ED  =  Emission  Delay  fer  the  secondary  station. 

ASF  =  Additional  Secondary  Phcise  Correction  for  point  P. 
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Equation  II-2  cjoiitdins  two  c»rrecticMis  vAiich  are  applied  to  Loran-C 
TD  readings.  The  velocity  of  electranagnetic  waves  in  air  under 

standard  cc»iditic»\5  is  changed  vhen  the  wave  propagates  over  a  surface. 

The  SeoCTTdcury  Phase  Oorrectian  Factor  (DTg  amd  DTj^  in  equatiCTi  II-2)  is 
a  time  correcticMi  for  the  change  in  velocity  caused  by  propagaticai  of 
electrcnvagnetic  waves  over  seawater.  This  change  in  travel  time  can  be 
calculated  by  the  following  errpirically  derived  equaticxis  (Ref.  2) : 

For  distances  greater  than  100  statute  miles  (T  >  537  nsec) : 

DT  =  Aq/T  +  Ai  +  A2*T  [II-3] 

vhere: 

EJr  =  Secondciry  Fhcise  Correction  for  all  seawater  path. 

Aq  =  129.04398  n-sec^ 

hi  =  -0.40758  /Ltsec 
A2  =  0.00064576438 

T  =  Travel  time  (Tg  or  in  equation  II-l)  in  /isec. 

For  distances  less  than  100  statute  miles  (T  <  537  /isec) : 

OT  =  Bq/T  +  Bi  +  B2*T  [II-4] 

vhere: 

Bq  =  2.7412979  /Lisec^ 

B]^  =  -0.011402  nsec 
B2  =  0.00032774624 

T  =  Travel  time  in  nsec. 
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In  equation  II-2  the  Additicxial  Secondary  Phase  Factor  (ASF) 
corrects  for  propagation  over  land  (actually  any  other  non-seawater) 
surfaces,  vhich  have  a  wide  variety  of  ocaiductivities  (see  Section 
III-C-2  or  Ref.  3  for  further  explanation  of  ASF  CorrecticxTs) .  The  ASF 
correcticxi  is  in  addition  to  the  Secondary  Phase  Correcticsi  for  all 
seawater  path. 

Secondary  Phase  Corrections  for  all  seawater  path  are  always  added 
in  before  calculation  of  the  geographic  position  from  TD's.  ASF 
corrections  are  niuch  snaller,  and  are  often  determined  from  tables 
which  require  a  dead  reckoning  positicxi  within  10'  of  latitude  and 
longitude  of  the  final  position.  If  no  CR  position  is  available,  the 
position  Ccilculated  from  TD's  corrected  only  for  all  seawater  path  is 
used  as  a  ER  position  to  determine  ASF  corrections.  When  used  in 
hyperbolic  mode,  the  ASF  correction  in  Equation  II-2  is  a  net 
correction  vhich  applies  to  the  particular  master-seccaidary  pair  (the 
difference  between  the  range-range  corrections  for  each  station) . 

Emission  Delay  (ED)  in  Equation  II-2  includes  a  coding  delay,  \(hich 
is  the  time  between  reception  of  the  master  pulse  at  the  secondary 
station  and  transmission  of  the  secondary's  pulse,  plus  one  way  travel 
time  from  the  master  to  the  secondary,  including  the  Secondary  Phase 
Correction  for  all  seawater  path  for  the  master-secondatty  baseline  path 
(Ref.  4). 

The  result  of  application  of  these  corrections  to  the  beisic  time 
difference  between  reception  of  master  and  secondary  signals  is  a 
reduction  of  the  time  difference  to  a  path  throui^  air  with  am  Index  of 
Infraction  of  1.000338  (Ref.4:p.4). 
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C.  Differential  CSarrocticns 

Differential  techniques  can  be  used  to  improve  the  absolute 
accuracy  of  loran-C  TD  readings.  Differentiad  Loran-C  uses  monitor 
receivers  located  at  fixed  points  in  or  neeur  the  aurea  vhere  inprovement 
is  needed.  The  most  coranon  method  is  to  use  TD  readings  of  the 
mcnitor  receivers  to  determine  the  magnitude  of  the  variations  in  the 
TD's  from  the  long  term  average  of  the  TD's  at  the  fixed  points.  These 
variations  are  then  applied  to  receiver  readings  in  the  area  of 
interest  before  geographic  positions  are  calculated  from  them. 
Corrections  are  xpdated  periodicailly  as  explained  later  in  this 
section.  (Ref .l:p.434) 

The  differential  correction  (including  geographic  correction)  to 
loran-C  readings  at  an  unknown  point  P2  with  a  monitor  receiver  at  a 
known  point  can  be  described  ky  the  following  equation: 

TD  at  P2  =  Received  TD  at  P2  +  Differential  Cor.  +  Geographic  Cor. 

or 

TD2C  =  TD2  +  (TDi  -  TDic)  +  [II-5] 

vhere 

TD2C  =  Corrected  TD  reading  at  Point  P2  (the  unJcxwn  point) . 

TD^c  =  long  term  average  TD  at  P^  (location  of  the  monitor) . 

TDj^  =  TD  reading  of  receiver  at  the  known  point  (Pj^) . 

TD2  =  TD  reading  of  the  receiver  at  the  unknown  point  (P2) . 

TC^  =  TD  variation  vhich  is  a  function  of  location. 


10 


For  practiced  purposes  it  is  possible  to  use  any  oorrectican  frtm 
real  time  to  a  correction  averaged  over  a  period  of  to  several  days. 
Ihe  method  used  in  this  experiment  (described  in  Secticai  IV-C)  was  to 
determine  a  regression  line  of  iKXiitor  receiver  ID  readings  at  the 
known  point.  A  computed  correction  was  applied  to  the  TD  reading  at 
the  unknown  point  based  on  the  time  of  the  TD  reading  at  the  unknown 
point  and  the  equaticxi  of  the  regression  line.  This  method  was  used 
because  CTily  one  mcaiitor  receiver  was  available  and  it  was  used  as  both 
a  monitor  at  the  known  point  and  at  the  unknown  points.  In  this 
experiment  a  correction  beised  on  reception  of  three  seoondcury  staticMis 
(rather  than  two)  was  applied  to  the  TD  readings  from  eaoh  secondciry 
station  before  the  regression  line  parameters  were  calculabed,  as 
explained  later  in  Section  IV-B. 

The  geographic  correction  -  that  part  of  the  difference  in  TD 
readings  at  points  and  P2  vAiich  is  a  due  to  the  differences  in  ASF 
corrections  at  the  two  points  -  is  a  pcurticulcu:  problem.  For 
geographic  points  at  sea,  tables  of  ASF  oorrecticsTS  within  10'  of 
latitude  and  longitude  are  available,  and  may  be  accurate  enough  to 
provide  corrections  within  the  necessary  accuracy.  Tables  of 
corrections  for  geograi*iic  points  on  land  are  not  published,  but  were 
obtained  for  this  eiperiment  from  the  Defense  Mapping  Agency 
Hydrographic/Topographic  Center. 

Test  results  of  differential  Loran-C  on  the  East  Coast  U.S.  chain 
have  shewn  accuracies  of  about  4.6  m  CEP  at  distances  up  to  140  km  from 
the  monitor  site  (Ref  .l:p.434) . 
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m.  (22IERAL  ANAI^SIS  OF  mCRS 


To  quantify  the  amount  of  error  and  compare  the  accuracy  of 
different  types  of  corrections  made  to  Icran-C  geograjiiic  position 
fixes,  it  is  necessary  to  provide  seme  basic  definitions.  The 
follcwing  sections  give  the  definitions  used  in  this  study  for  the 
types  of  accuracy  and  measures  of  accuracy.  One  and  two  dimensional 
sources  of  error  in  loran-C  positioning  are  listed  and  approximate 
values  for  the  error  cire  given  where  possible.  Error  reduction  based 
on  differential  correction  is  described. 

A.  Types  of  Accuracy 

There  are  several  ways  to  specify  the  accuracy  of  electronic 
navigation  systems.  The  following  definitions  of  three  types  of 
accuracy  have  been  found  useful  in  this  study  (from  Ref.  5) ; 

Absolute  Accuracy;  For  purposes  of  this  study,  absolute  accuracy 
is  a  mecisure  of  the  ability  to  determine  geographic  position  (in  terms 
of  latitude  and  longitude  or  other  fixed  coordinate  system)  fron  a 
given  set  of  Loran-C  coordinates  (Ref.5:p.l72) .  The  error  includes  all 
types  of  systematic  eind  randem  error  associated  with  the  system. 

Absolute  accuracy  in  this  experiment  refers  to  the  accuracy  of  the 
calculation  of  the  positions  of  the  "unknewn"  points  based  on  receiver 
TD  readings,  three  station  corrections,  and  differential  corrections, 
ccnpared  to  the  knewn  published  positions  of  those  points. 
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Repp.al-.ahle  Aocuracv;  Ihe  ccmtiDn  navigation  definition  of 
repeatable  accuracy  is  the  ability  to  return  to  the  same  point  using  a 
given  set  of  electronic  coordinates  (Ref  .5:p.l72) .  Error  is  measured 
in  linear  units.  This  differs  from  absolute  accuracy  in  that  the 
position  itself  need  not  be  known  in  terms  of  seme  other  coordinate 
system.  Repeatable  accuracy  is  normally  better  than  absolute  accuracy 
because  the  part  of  systematic  error  vAiich  is  a  furction  of  position 
has  been  removed  (the  position  remains  the  same) ,  leaving  only  rcmdcm 
error  and  time  dependant  systematic  error.  Repeatable  accuracy  in  the 
context  of  this  eiqjeriment  is  a  measure  of  the  precision  of  measurement 
of  the  TD's,  vhich  is  most  affected  by  the  electronic  circuitry  of  the 
transmitters  and  receivers  involved  (over  short  periods  of  time  in 
which  terrain  and  atmospheric  factors  do  not  change  significantly) . 

Repeatable  accuracy  is  similar  to  the  accuracy  eigjected  fron 
differential  corrections.  The  difference  between  the  two  is  that 
since  repeatable  accuracy  refers  to  the  ability  to  return  to  the  same 
point  using  the  same  set  of  received  TD  readings,  there  is  no 
difference  in  the  ASF  (terrain)  correction  but  there  may  be  some 
difference  due  to  V2triability  of  TD's  with  time. 

Differential  corrections  have  different  ASF  corrections 
at  different  points.  These  must  be  calculated  and  applied 
to  the  TD  readings.  At  a  single  point  the  ASF  corrections 
will  probably  not  change  over  the  time  period  of  interest, 
but  over  longer  periods  some  error  will  be  introduced.  The 
uncertainty  in  ASF  corrections  could  be  a  significant 
source  of  error. 
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Relational  Anairacy:  The  accajracy  with  which  it  is  possible  to 
determine  the  position  of  one  observer  in  relaticxi  to  another  observer 
vising  the  same  system  (Ref.5:p.l72) .  A  form  of  relaticMVil  aocuracy 
enters  into  differential  oorrecticais  when  tvro  different  receivers  are 
vised.  Even  at  the  same  point  at  the  same  time  there  will  usucilly  be  a 
sli^t  difference  in  the  loran-C  readings  of  different  receivers.  This 
is  caused  by  errors  aussociated  with  the  receiver  circuitry  cuid  randcm 
errors  over  the  sli^t  differences  in  the  time  at  v*iich  the  readings 
are  made  by  each  receiver,  since  all  other  factors  are  the  same.  The 
same  type  of  error  exists  with  two  different  receivers  s^erated  by 
seme  distance  (ais  with  differentiad  corrections) ,  adthou^  this  type  of 
error  is  usuadly  small  cenpared  to  other  sources  of  error. 

B.  Meaisures  of  Aocuracy 

In  addition  to  the  types  of  accuracy  listed  above,  it  is  necessary 
to  specify  hew  aocuracy  is  to  be  meaisured.  This  cam  be  viewed  ais  am 
estimate  of  the  aonount  of  error.  The  following  eune  definitiems  of  the 
basic  methods  used  here  of  specifying  the  measures  of  accuracy  of 
position.  Their  vasefvdness  and  limitations  in  this  application  eire 
also  described. 

Root  Mean  Square  fPM5)  Error;  This  is  the  linear  (one  dimensicmal) 
error  in  the  location  of  an  IDP  (Ref.5:p.l72) .  It  can  be  specified  in 
lineau:  units  or  time  units  (meters  amd  microseconds  are  used  here) . 
The  normal  specification  of  RMS  error  is  the  number  of  meters  or 
microseconds  r^resenting  sane  vhole  number  of  standaund  deviations  from 
the  mean.  In  the  application  used  here,  the  variation  in  TD  readings 
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frcm  each  master-secondary  pair  can  be  ejqsressed  as  a  steuidard 
deviation  in  microseconds,  vhich  corresponds  to  a  displacement  of  the 
LDP  by  a  corresponding  nuitber  of  meters  measured  normal  to  that  line  of 
position. 

Circulcir  Error  Prcfcable  (CEP^ ;  This  is  the  two  dimensional  error 
vhich  specifies  the  radius  of  a  circle  within  which  there  is  a  50% 
probability  of  being  located,  based  on  positions  calculated  fran  the  TD 
readings  of  the  receiver.  The  estimated  position  of  the  point  (based 
on  latitude  eind  longitude  or  other  coordinate  systems)  and  the  radius 
of  the  circle  are  used  to  specify  this  measure  of  accuracy.  Circles 
other  than  the  standard  50%  circle  can  be  used  (67%  eind  90%  eure 
common).  (Ref.5;p.l73)  The  main  advantage  of  using  CEP  (rather  than 
the  error  ellipse  described  below)  is  that  direct  comparison  between 
accuracies  of  different  systems,  techniques,  or  at  different  points  can 
be  made.  CEP  in  this  experiment  is  used  as  a  measure  of  the  precision 
of  the  TD's  indicated  by  the  receiver  at  a  fixed  point  eind  the 
precision  of  corrected  TD's  based  on  these  uncorrected  TD's. 
Corrections  for  the  absolute  position  (the  bias  of  the  TD's  from  true 
readings)  is  of  secondary  inportance  in  this  application. 

Circular  error  estimates  in  this  experiment  are  based 
on  the  radius  which  50%  of  the  data  are  greater  and  50%  are 
less,  as  calculated  by  a  sorting  program. 
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Error  Ellipses;  Eiror  ellipses  are  the  preferred  way  of 

describing  the  accuracy  of  electrmic  navigation  systems  for  many 
purposes.  As  with  circular  error,  accuracy  is  ^jecified  in  terms  of 
the  probability  of  an  estimate  of  geographic  positicxi  feilling  within  a 
gecmetric  figure  centered  on  the  actual  geographic  position  (ein  ellipse 
rather  thari  a  circle  is  used)  (Ref.  6) .  Unlike  the  circle  used  in  OEP, 
cm  ellipse  provides  additicHial  informaticMi  about  the  directiaieLl 
characteristics  of  the  accuracy,  vhich  can  be  very  inportant  in  scroe 
applications.  The  specifications  for  elliptical  accuracy  are  usually 
the  lengths  of  the  major  cind  minor  axes  and  the  orientaticai  of  the 
ellipse  (such  as  the  azimuth  of  the  direction  of  the  major  axis) .  A 
disadvantage  of  using  error  ellipses  in  this  experiment  is  t2ie 
difficulty  of  directly  oenpetring  the  accuracy  of  different  ellipses 
with  different  orientations.  For  that  reason  most  of  the  accuracy 
oenparisons  in  this  study  have  been  specified  in  terms  of  circular 
error.  Error  ellipses  have  been  used  vhere  appropriate. 

Error  ellipses  are  related  to  the  orientation  of  two  intersecting 
LOP'S  by  the  following  equations  (Ref. 6:p. 72-74) : 


tan  20  = 


sin  20 


oi^(coe  20)  +  02^ 


[III-l] 


2  _ 


2  sin^0 


[III-2] 
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[III-3] 


2  sin^^ 


where: 

0  =  the  ^nqle  of  the  major  axis  of  the  ellipse 

from  the  direction  of  the  LOP  with  smaller 
variance,  toward  the  LOP  with  greater  variance, 
in  the  direction  of  the  smaller  angle. 

=  variance  of  the  IDP  with  smaller  variance. 

02^  =  Vciricince  of  the  IDP  with  the  greater  variance. 

0  =  the  smaller  angle  between  LOP ' s . 

Ihese  equations  were  used  to  calculate  the  dimensions  and 
orientation  of  the  error  ellipse  for  Range-7,  the  ''kncwn”  point  used 
for  differential  corrections  (see  Figure  V-3) . 


C.  Souroes  of  Linear  Error  in  loran-C  Time  Differences 

The  position  fixing  accuracy  of  the  LDRAN-C  system  depends  on  a 
number  of  factors.  Ihese  can  be  divided  into  two  categories 
(Ref  .7:p.l8) :  tenporal  errors  affecting  the  repeatability  of  the 

system  (errors  which  vary  with  time) ,  cind  systematic  (time  inveiriant) 
errors.  Both  of  these  types  of  error  may  also  be  functions  of 
geographic  position. 

The  effect  of  these  errors  is  to  produce  TD  readings  at  the 
receiver  which  differ  from  those  which  would  be  expected  on  a  perfect 


hypertxjlic  lattice  cai  the  surface  of  a  perfect  ellipsoid,  resulting  in 
error  in  the  calculation  of  the  geographic  positiOTi.  A  lineau:  error  of 
1.0  Msec  in  time  difference  is  approximately  equivcilent  to  150  m  on  the 
baseline,  with  caie  dimensicxieLL  error  at  other  points  being  i?)  to  abcxit 
ten  times  greater,  d^jending  on  the  locaticxi  within  the  gecmetry  of  the 
system. 

The  main  sources  of  error  in  the  determination  of  TD's  used  to 
ccilculate  geogr^ihic  position  can  be  classified  as  follorfs  (Ref. 8); 

•  Transmitter  synchronization  timing  errors. 

•  Errors  associated  with  propagation  effects. 

•  Atmospheric  noise. 

•  Receiver  errors. 

1.  Transmitter  Synchronization  Timing  Errors 

Since  time  differences  between  reo^>tion  of  a  itaster  pulse  and 
secondary  pulses  are  used  to  calculate  geographic  position  in 
hyperbolic  loran-C,  it  is  essenticil  that  the  synchronization  of  timing 
between  transmissions  from  the  master  and  secondaries  be  hi^ily 
precise. 

In  order  to  acocnplish  this,  each  station  is  individually  timed 
using  cesium  clocks  vhich  drive  pulse  and  grotp  r^etition  interval 
timing  circuits.  Ind^jendent  monitor  receivers  are  used  to  provide 
data  to  ccnputers  vhich  adjust  the  relative  timing  between  master  and 
secondaries  of  the  chain  in  small  st^s  (usually  10  to  20  ns) . 
Absolute  timing  of  the  master  stations  is  adjusted  to  c^jerate  in 
coordination  with  Universal  Time,  Coordinated,  (Ref .8:p.  1129) 
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Estimates  of  timing  error  due  to  transmitter  synchronization 
are  on  the  order  of  0.05  /isec  (Ref.9:p.537) .  Actual  data  from  eui 
operaticaieil  Loran-C  chain  showed  time  differences  controlled  within 
0.024  jusec  of  the  stendaud  vailue,  with  a  standard  deviatic»i  (in 
ed^solute  time)  of  0.032  nsec  (Ref.  10' p. 227) .  A  transmitter 
synchrmization  error  of  0.01  ^sec  is  believed  to  be  attainable  if 
ground  mcxvitor  data  is  used  to  remove  the  meein  synchrOTiizaticxi  error 
(Ref.ll:p.401) . 

2.  E^rrors  Associated  %n.th  Propagation  Effects 

This  includes  most  types  of  error  external  to  the  electrcxiics 
of  the  system,  and  is  the  main  source  of  error  in  the  determinaticsi  of 
geographic  position. 

Factors  vhich  influence  groundwave  propagation  include  the 
following  (Ref .12; pp. 39-53  and  Ref. 13 ;pp. 173-187) : 

•  Electrical  Conductivity  of  the  surface  material 

•  Dielectric  Constant  of  the  surface 

•  Index  of  Refraction  of  the  air 

•  Lapse  Rate  of  the  Index  of  Refraction  with  altitude 

These  factors  are  ocnplex  functions  of  several  variables,  the 
main  ones  being: 

•  type  and  conposition  of  the  surface  over  which  the  wave  travels 

•  topography  of  the  surface 

•  meteorological  conditions 


At  any  particular  geographic  positiCTi  the  effect  of  these 
variables  may  or  may  not  be  considered  a  function  of  time,  d^5ending  cMi 
the  time  periods  of  interest. 

Also  of  importance  is  the  effect  of  these  factors  on  ID  as  a 
function  of  geographic  position.  For  each  different  geogre^iuc 
position,  the  electromagnetic  waves  from  each  transmitter  must  traverse 
a  different  route,  with  different  types  of  surface,  tcpography,  and 
atmospheric  oonditicMTS,  to  reach  the  receiver. 

The  effects  of  irregular  terrain  are  practically  time 
invariant,  while  ground  conductivity  and  dielectric  constant  may  vary 
seasoned ly  (due  to  cheinges  in  moisture  content  of  the  soil,  changes  in 
vegetation,  or  other  factors) .  CSaserved  cyclic  annual  variations  cai 
the  order  of  0.5  /isec  aune  probably  due  to  seasonal  Vcuriations  in  ground 
conductivity  eind  dielectric  constant  and  to  seasonally  vaiying 
atmo^heric  ocHiditicaTS . 

Meteorological  conditions  sufficient  to  cause  noticeable 
variations  in  TD  readings  may  change  in  very  short  time  periods  (days 
or  hours) .  Short  term  vzariaticais  due  to  meteorological  effects  (cold 
front  esri  warm  front  passage  through  the  line  of  propagation)  on  the 
order  of  0.025  Msec  have  been  observed.  These  correlate  with  chcinges 
in  the  dry  adieibatic  Ic^sse  rate  of  the  index  of  refraction. 
(Ref. 14 :p. 1119) . 

For  most  differential  loran-C  purposes  it  is  necesscury  to 
consider  cnly  those  tenporal  propagation  effects  vhich  occur  in  short 
time  periods  (hours  or  days  at  most,  for  most  applications) .  The 
variation  of  propagation  effects  with  geographic  position  must  also  be 
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considered  for  precise  positioning.  Many  of  these  variaticxis  can  be 
corrected  by  a  corputation  of  the  AdditicMicLL  Seccxidary  Phase  Factor 
(ASF). 

The  Additional  Seocxidary  Phase  Factor  (ASF)  is  a  measure  of  the 
systematic  error  resulting  from  differences  in  velocity  of  electrcmag- 
netic  waves  over  ground  with  different  ocxiductivities  and  dielectric 
constants.  The  Defense  Mapping  agency  uses  a  ccoputerized  system  of 
maps  with  estimates  of  ground  ccnductivity  (effectively  including 
dielectric  constcint)  over  large  areas,  and  ceilculates  ASF  corrections 
in  the  loran-C  coverage  areas.  These  corrections  are  published  as 
tables  of  ASF  corrections  to  be  used  in  correcting  Loran-C  TD  recKiings 
from  receivers. 

ASF  corrections  can  be  as  leirge  cis  4  microseconds  (Ptef.3). 
The  corrections  provided  in  the  Defense  Mapping  Agency  tables  require 
estimates  of  the  geographic  locaticai  of  the  receiver  to  the  nearest  10* 
of  latitude  and  longitude.  ASF  correctiOTs  based  on  the  known 
positions  of  the  survey  points  used  in  this  ejqjeriment  were  provided  by 
EMA  using  a  ccnputer  program  vhich  calculates  the  corrections  using  the 
Millington  Method  (Ptef.3;p.V) . 

3 .  Atmospheric  Noise 

The  effects  of  atmospheric  electrical  noise  are  reduction  of 
the  range  at  vhich  Loran-C  signals  can  be  effectively  received  and  an 
increase  the  error  in  TD's  ccnputed  by  receivers  in  marginal  recepjtion 
areas.  This  is  primarily  d^jendent  on  the  circuitry  of  the  receiver. 
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In  roarginea  areas  the  effect  of  atmospheric  noise  will  be  an 
increeise  in  the  stcindcirtl  deviation  of  the  ID's  and  increased  difficailty 
in  locking  cxi  and  maintaining  recepticxi  of  the  weaker  signals. 
Atmo^heric  noise  varies  with  time  eind  geographic  locatiCTi,  with  a 
significant  diumed  variaticn  (noise  levels  are  higher  at  ni^t) .  In 
this  experiment,  significant  effects  due  to  local  sources  of  electrical 
noise  were  noticed  in  the  Mcxiterey  h2Ud»r  eunea. 

4.  Receiver  Errors 

Receiver  errors  eire  of  two  main  types: 

•  timing  error  caused  by  limitations  of  the  intemcil  clock 

circuitry  of  the  receiver 

•  errors  caused  by  processing  time  of  the  receiver 

For  a  receiver  at  a  fixed  geographic  location,  processing  time 
is  not  a  source  of  error.  Ihe  receiver  will  be  at  the  same  location 
regardless  of  the  time  required  to  output  the  TD's  or  geographic 
coordinates.  The  main  source  of  receiver  error  at  a  fixed  locaticxi 
(other  than  the  noise  considerations  discussed  above)  is  the  timing 
error  in  the  internal  clock  circuitry.  Loran-C  receivers  vary  greatly 
in  their  timing  errors.  Ihe  Racal  Megapulse  Monitor  receiver  used  in 
this  experiment  has  a  specified  timing  error  of  0.02  jusec  or  better. 

For  mobile  receivers  the  processing  time  of  the  receiver 
circuitry  can  be  important.  After  the  signal  is  received  it  must  be 
processed  and  either  stored  in  the  form  of  electronic  information  or 
visually  displayed.  Error  is  introduced  if  the  receiver  moves  during 
this  time  interval. 
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Ctoviously,  the  anount  of  error  in  the  determinaticai  of 
geographic  position  from  this  source  depends  cxi  the  processing  time  and 
the  speed  of  the  receiver.  Receivers  vAiich  output  data  in  the  form  of 
geogrcphic  coordinates  nust  calculate  latitude  and  Icaigitude  from  the 
TD  readings  of  the  receiver,  requiring  more  processing  time  and 
resulting  in  greater  error  from  this  source. 

At  the  speeds  normally  used  in  hydrographic  surveying 
(typically  less  than  20  kt)  the  processing  time  is  insignificant  for 
practically  all  receivers.  For  data  vthich  is  stored  in  the  form  of 
receiver  output  TD's  for  later  correction  and  conversion,  as  with  the 
Racal  Megapulse  receiver  used  here,  error  should  be  insignificant  even 
at  much  hi^er  speeds. 

D.  Geometric  Dilution  of  Ptecisicn 

This  source  of  error  in  geographic  position  results  from  the 
crossing  angle  of  the  loran-C  IDP's  (lines  of  position) .  Figure  IV-1 
shews  two  hyperbolic  IDP's  with  standard  deviaticais  and  02  (in  time 
difference)  and  crossing  eingle  a.  The  optimum  crossing  angle  of  any 
two  IDP's  is  90°,  resulting  in  minimum  error  in  geographic  position  for 
a  given  set  of  standard  deviations. 
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E.  Summatic^  of  Errors 


1.  Absolute  Aocuracy  vs.  Frecisixxi 

The  absolute  acxairacy  of  an  electrcxiic  positicxiing  system  can 
be  viewed  as  a  measure  of  its  ability  to  give  the  geogre^iiic  position 
of  the  receiver.  Precisicai  is  a  measure  of  the  variation  of  individual 
time  difference  readings.  If  the  receiver  is  set  vp  at  a  point  for  a 
sufficient  amount  of  time,  many  of  the  randan  errors  can  be  averaged 
out;  the  average  of  a  large  nurrber  of  readings  can  provide  a  basis  for 
determining  absolute  aocuracy. 

In  practice  this  is  more  difficult  than  it  would  seem,  since 
sc3me  factors  vJiich  do  not  vary  significantly  over  short  periods  of  time 
(hours  or  days)  may  vary  significantly  over  loiger  periods  (weeks  or 
months) .  It  is  necessary  to  consider  the  time  periods  involved  in 
order  to  determine  vAiich  factors  to  consider  random  and  ViAiich  to 
consider  "bias”  in  the  system. 

For  this  experiment,  the  sources  of  random  error  are  coisidered 
to  be  transmitter  timing  errors,  some  iteteorological  oanditioTs, 
atmospheric  noise,  and  most  receiver  errors.  Errors  caused  by  long 
term  "bias"  of  the  system  are  due  primarily  to  ground  conductivity  and 
type  of  surface,  sane  meteorological  effects,  and  any  non-random  bieis 
in  the  timing  of  the  transmitters. 

Taking  these  factors  into  consideration,  data  tciJcen  at  the 
"kncwn"  point  (Range-7)  of  the  differenticil  system  was  processed  to 
remove  random  error  by  averaging  over  a  period  of  time,  vrtiile  long  term 
error  was  removed  by  caiparison  of  the  TD's  to  their  expected  values 
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cSerived  from  calculaticaTS  provided  by  oorputer  cxirputaticxis  of  the  TD 
at  that  point. 

Errors  which  may  change  within  the  time  frame  of  the  experiment 
were  reduced  by  a  linear  regression  of  the  trend  of  data  taken  at  the 
known  point  and  projected  to  the  correction  at  the  unknown  point 
mathematiCcilly. 
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IV.  IHREE  STATION  AND  DIFFEPENITAL  OC^RECnOIS 


A.  Basic  Gecmetxy  of  loran-C  Hypertx>leis 

Hus  secticxi  provides  a  sumrury  of  the  information  neoesseiry  for 
the  analysis  of  Loran-C  oorrecticxTS  for  three  seocaidary  staticsis  and 
for  differential  Loran-C  oorrecticMTS  presented  in  p)arts  IV. B  cind  IV.  C, 
adapted  fron  Laurilla  (Ref.l)  emd  others  as  noted. 

1.  Lanewidth 

Hie  baseline  of  a  loran-C  roaster-secondary  pair  is  the  shortest 
path  over  which  the  transmitted  signals  may  travel  between  the  roaster 
station  and  the  secondary  station.  On  a  spherical  earth  model  it  is 
the  eox:  of  the  great  circle  through  the  positions  of  the  roaster  cuid 
secxindary  stations.  On  eui  ellipsoid  it  is  the  geodesic  from  the 
master  to  the  secondary  station.  In  the  set  of  hyperbolas  of  time 
differencje  for  the  master-secondeiry  pair,  the  s^aration  between 
hypertxDlas  is  a  minimum  constant  value  along  the  baseline. 

A  lane  is  the  distance  between  two  adjacent  hyperbolas  of  unit 
time  difference.  Hie  unit  normally  used  is  the  /isec  (10“^  sec).  In  a 
series  of  hyperboleis,  as  shewn  in  Figure  II-l,  it  is  obvicxis  that  the 
lemewidth  varies  with  location.  At  emy  point,  lanewidth  can  be 
viewed  eis  the  inverse  of  the  gradient  of  the  ratio  of  the  change  in 
time  difference  to  the  change  in  disteince  (Ref.  15),  or: 

L  =  {|(dTD/dl) I )~^  [IV-1] 
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MOST  FROBABIE  POINT 


Figure  IV-1.  Typical  intersection  of  two  lines  of 
position.  Most  probable  point  is  at  the  Intersection  of 
the  two  LOP'S. 


Figure  IV-2.  Typical  intersection  of  three  lines  of 
position.  When  three  or  more  secondary  stations  can  be 
received  at  a  geographic  position,  calculated  position 
depends  on  the  LOP ' s  which  are  used  to  make  the 


The  lanewidth  at  any  point  in  the  set  of  hypejtolas  can  be 
expressed  in  terms  of  the  lane  width  on  the  kaseline.  Beiseline  lane 
width  can  be  eiqiressed  by  the  follcwing  equation: 

Ib  =  -^  [IV-2] 

vAiere: 

Ljj  =  Lanewidth  on  the  baseline  (nebers/nsec) 

c  =  Speed  of  electrcroagnetic  waves  on  the  baseline  (m/i.Lsec) 
(note  that  time  units  are  in  the  unit  lanewidth) 

Lanewidth  at  any  other  point  (not  on  the  baseline)  can  be 
expressed  by  equation  IV-3a,  where  p  is  the  angle  subtended  by  the 
baseline  from  the  point.  Note  that  this  angle  is  the  absolute 
difference  between  the  azimuths  frcm  the  point  to  the  roaster  and  fran 
the  point  to  the  secondary  station. 

L  =  lte{sin(;0/2))-l  [IV-3a] 

or 

L  =  IfcG 

where: 

L  =  Lanewidth  at  a  point  (per  unit  of  time) 

=  Lanewidth  chi  the  baseline  (per  unit  of  time) 

P  =  Angle  subtended  by  the  baseline  from  the  point 
G  =  The  Lane  Expansion  Factor,  G  =  {sin(/3/2)  }"^ 

Note  that  the  lane  expansion  factor  depends  only  on  the  angular 
separation  of  the  master  and  the  secondary  stations.  This  ctiviously 
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d^iends  an  the  d-isd-aiw*  of  the  point  from  the  baseline  center  and  the 
angular  separation  of  the  baseline  zuid  the  vector  from  the  baseline 
center  to  the  point. 

Equation  IV-3a  provides  a  mecins  of  calcul  ating  the  normal 
distance  between  emy  two  hyperboleis,  provided  th^  are  near  encu^  that 
their  lanewidths  are  not  significantly  different: 

d  =  Lt  [IV-3b] 


vAiere: 

d  =  distance  in  meters  between  the  two  hyperboleis. 

L  =  lanewidth  at  the  point  of  interest  (in  ra/^sec) . 

t  =  time  difference  between  the  two  hyperbolas  in  jusec. 

2.  Direction  of  Hyperbola 

The  direction  of  a  hyperbola  at  any  point  on  the  hyperbola  is 
cilong  the  teingent  at  that  point,  directed  away  from  the  baseline  of  the 
master  and  secondary  stations  used  to  define  the  hyperbola.  For 
applications  in  this  ejqjeriment,  Laurilla  (Ref  1)  provides  the  follow¬ 
ing  useful  definition  of  the  direction  of  a  hyperbola:  "...  the 

direction  of  a  hyperbola  at  any  point  coincides  with  the  bisector  of 
the  angle  formed  by  the  lines  joining  the  point  to  the  pair  of 
stations."  (Ref.l:p.92) 

For  purposes  of  this  experiment,  it  is  convenient  to  ei^ress 
the  direction  of  the  hyperbola  by  its  azimuth.  This  azimuth  can  be 
expressed  in  terms  of  the  azimuths  (from  north)  of  the  master  and 
secondary  stations  by: 
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“h  =  +  “s)/2  (^-^lere  1 0^-012 1  >180*)  [IV-4a] 

Oh  =  180*  +  (£%  +  o^)/2  (where  |aiii-as|<180')  [IV-4b] 

vtiere: 

ojj  =  aziinuth  of  the  direction  of  the  hypertxDla. 

Ojj,  and  OEg  =  eizinuths  of  the  master  and  secaidary  from  the  point. 


3.  i^plicahicns  of  Gecmetry  of  Hyperbolas  to  E^rror  Aralysis 

The  geometric  description  of  loran-C  hyperbolas  in  parts  IV. A.  1 
and  rv.A.2  is  particulzarly  inportant  in  the  cincdysis  of  error  eind  the 
derivation  of  corrections  used  in  this  experiment.  Seme  of  the  direct 
applications  are  outlined  in  the  following  paragraphs. 

The  error  ejpressed  as  an  error  ellipse  requires  the  directicans 
of  hyperbolcis  in  order  to  determine  the  orientation  of  the  ellipse. 
The  dimensions  of  the  ellipse  require  the  conversion  of  variances  of 
time  differences  of  each  hyperbola  to  variances  of  distance  in  the 
major  euid  minor  axis  of  the  ellipse.  Standard  deviations  of  time 
differences  can  be  converted  to  distances  by  the  following  equations: 

0(3  =  [IV-5a] 
Od^  =  [IV-5b] 
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The  cx)rrecticMi  for  three  seccxidary  stations  uses  variances  in 
time  differences  converted  to  distances  as  an  input,  using  equation 
IV-5b.  Converting  time  difference  to  distance  in  this  way  correctly 
wei^ts  the  precision  of  the  lOP's  vAien  determining  the  correction 
factor  based  on  gecnetric  ccH^ideraticHis. 

B.  CSorrecticn  for  Three  Seocndary  Stations 

In  situaticMTs  where  three  seocxxiary  staticais  can  be  received,  the 
result  is  three  LDP’s  vhich  do  not,  in  general,  intersect  at  one  point 
(see  Figure  IV-2)  viien  they  are  used  to  ceilculate  the  geographic 
position,  althou^  they  must  in  fact  intersect  at  the  common  point  at 
vhich  the  readings  were  taken  concurrently.  Uncier  usuad  navigation 
conditions  the  two  best  LDP's  are  used  to  determdne  a  geographic  fix. 
A  dra\^ck  of  this  method  is  that  it  fails  to  consider  all  of  the 
available  information  (i.e.  one  of  the  available  lines  of  position  is 
not  used  in  the  calculation  of  the  geographic  position) .  By  utilizing 
all  three  LDP's  it  should  be  possible  to  improve  the  accuracy  of  the 
geographic  fix. 

Sections  B.  1  and  B.2  will  develop  a  mathematical  adjustment  which 
can  be  used  to  calculate  the  most  probedjle  cxjrrected  TD  readings  of  a 
Loran-C  receiver  used  as  the  fixed  receiver  in  a  differential  Loran-C 
system  utilizing  three  LDP's  (i.e.  receiving  the  master  and  three 
secondary  stations) .  When  the  oorrecrt:ion  is  applied,  the  three 
corrected  TD's  will  intersect  at  the  most  probable  point,  based  on  long 
term  average  values  of  the  TD's.  The  derivation  of  this  adjustment 
cxjuld  take  emy  of  several  possible  forms,  but  a  primary  consideration 
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here  is  the  need  for  a  form  which  is  eeisily  e^^plicable  to  oonpiter 
processing.  Secticxi  B.3  describes  a  method  of  applying  the  corrections 
for  ncMT-differenticil  use,  using  only  one  receiver  and  with  no  kncwn 
fixed  point. 

Section  B.4  will  develc^  regressicxi  equaticxis  for  ceilcalating  the 
corrected  TD's  directly  from  time  for  use  as  a  differential,  correction 
for  a  reanote  receiver.  Note  that  the  methods  developed  in  these 
secticais  aure  primeurily  for  use  in  the  differential  loran-C  conputations 
necessauY  for  this  esqaeriment. 

1.  Geometric  CScnsideraticns  Used  to  Develop  a  Correction  Factor 

The  following  assuitptions  are  made  for  this  derivation: 

•  Each  IDP  has  a  unique  direction  (described  in  Section  IV. A. 2) 

•  Long  term  average  TD's  (designated  TDl^y,  TD2av»  arid  TBS^v)  fof 

for  the  known  point  are  available 

•  Vatriances  of  the  TD  vailues  (designated  ,  a-^ ,  and  are 

known  or  can  be  calculated 

•  Variations  of  TD's  from  their  average  values  is  in  the  form  of  a 

Gaussian  distribution 

•  Over  the  relatively  small  area  being  considered,  the  surface  can 

be  considered  a  plane  and  IDP's  for  the  same  master-secondary 

pair  can  be  considered  parallel  strai^t  lines 

Several  things  should  be  noted  about  these  assunptions.  Long 
term  averaging  times  for  this  ejqjeriment  will  be  on  the  order  of  a  few 
hours  during  acquisition  of  data;  for  purposes  of  the  derivation, 
standcud  deviations  cure  converted  between  time  units  (nsec)  and 
distance  units  (m)  using  Equation  IV-5a;  and,  the  assunptions  of  a 
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Gaussian  distribution,  plane  surface,  and  pcu:3llel  LDP's  eire  reasonable 
for  this  type  of  data  (Ref.6:p.8).  When  two  LDP's  with  different 
standard  deviations  intersect  (see  Figure  IV-l) ,  the  point  of 
intersection  is  the  most  probable  geographic  position  (Ref ,6:p. 8-12) . 
Ihe  probability  of  the  location  of  the  true  position  is  best  described 
using  ein  error  ellipse  as  described  in  Section  III-B. 

A  typical  situation  is  shown  in  Figure  3,  a  set  of  three 
concurrent  LDP's  intersecting  at  three  different  points.  Point  P  is 
the  location  of  the  monitor  receiver  used  to  determine  differential 
corrections,  with  TDl^y,  TD2ay,  and  being  the  long  term  average 

TD  values  at  that  kncwn  point  (note  that  in  figure  IV-3,  the  LDP's  of 
TDl^vi  and  are  not  shewn,  but  are  parallel  to  TD^^,  TD2, 

and  TD3,  respectively,  and  pass  through  point  P,  as  shewn  in  Figure 
IV-4) .  Althou^  the  long  term  average  LDP  values  intersect  at  point 
P,  tricingle  P12P13P23  repi^'esents  a  triangle  of  fixes  formed  by  the  set 
of  TD  values  vhich  may  occur  from  ciny  concurrent  set  of  TD's,  vhich  are 
unlikely  to  equal  the  long  term  average.  Each  set  of  TD  readings 
(which  occurs  in  the  data  every  few  seconds)  will  result  in  a  new 
triangle  of  fixes. 

The  vectors  r^,  r2,  and  r3  represent  the  difference  between  the 
individual  TD  values  in  a  concurrent  set  of  data  and  their  respective 
long  term  TD  values.  The  r  vectors  are  at  ri<^t  angles  to  both  the 
long  term  and  current  LDP's.  For  exanple,  if  the  long  term  average  of 
TDj^  (TDlgy.)  is  29939.200  and  the  value  of  TD^^  in  a  set  of  TD  readings 
currently  being  considered  is  29939.300,  the  value  of  r^  for  the 
current  set  of  TD  readings  is  0.100  (all  in  ij,sec)  .  When  converted  to 
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meters  using  equation  IV-3b,  this  is  the  distance  between  the  average 
TD  values  (LDPigy)  and  current  TD  values  (TOi) .  These  LDP's  can  be 
considered  cis  parallel  straight  lines  over  the  relatively  small 
distcinces  considered  for  this  correction.  Note  that  there  are  three 
LDP's  in  each  set  of  concurrent  data  representing  TD  vcilues  for  three 
master-secondary  pairs.  Ihe  mcister  station  is  ccrmon  to  each  of  these 
pairs. 

Within  the  triangle  of  fixes  P12P13P23  by  the  three 

LDP's  of  a  set  of  data,  point  P*  represents  the  most  prt±iable  point  of 
the  true  geographic  position  indicated  by  that  particular  concurrent 
set  of  data  corpeired  to  the  long  term  average  data.  The  distances  of 
P'  fron  each  LDP  (the  lengths  r^')  are  related  to  the  relative 
precision  of  each  LDP,  P'  being  closer  to  LDP's  with  hi^er  precision 
(to  be  discussed  in  more  detail  in  the  next  section) . 

The  problem  must  be  solved  in  general  terms  so  that  any  set  of 
TD's  in  the  data  can  be  corrected  to  a  new  set  of  TD's  which  intersect 
at  the  most  probable  point.  Since  the  triangles  formed  by  all  sets  of 
data  are  similar  triangles,  the  lengths  of  r^',  r2',  and  r3',  will  have 
the  same  proportion  for  each  set.  Ihe  corresponding  r^'  values  between 
sets  of  data  will  be  related  by  a  ccsranon  factor,  which  will  be 
designated  as  the  Correction  Factor,  k. 

Other  information  which  should  be  noted  about  Figure  IV-3 
include  the  following: 

•  The  following  sets  of  lines  can  be  considered  parallel  over  the 
area  represented  by  the  figure:  TDl-v  II  TD^  ;  TD2av  II  TD2 ; 

TD3av  i|  TD3;  ri  II  ri'  ;  r2  f  r2'  ;  and  r3  |1  r3'. 

(continued  on  next  page) 
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•  'Ihe  following  sets  of  lines  are  perpendicular:  and  J- 

TDlg^  and  ;  r2  and  r2'  -L  TD2av  TD2  ;  r3  and  r3'  J. 

and  TD3. 

•  Ihe  azimuths  (from  north)  of  the  directions  of  the  LDP's 

represented  by  TDl^y^,  TDZg^yf,  and  cure  a^,  02,  eind  03, 

respectively.  Ihese  are  tiie  same  as  tdie  azimuths  of  TD^,  TD2, 
and  TD3,  respectively. 

•  Ihe  following  angles  of  triangle  P12P13P23  defined  by: 

“12  “  “2  ~  “1  •  ~  (^3  ~  ^1)  i  ®23  ~  “3  “  “1* 

These  angles  will  be  the  same  in  the  simileir  trieingles  formed 
by  any  lengths  of  ri,  r2,  and  r3. 


•  LDP's  1  and  3  have  the  largest  absolute  angle  (+  or  -)  between 
them,  and  LDP  2  is  the  intermediate  angle. 


The  next  section  will  prtjvide  a  mathematical  derivation  of  the 
e>pression  for  a  correction  factor  to  determine  the  lengths  of  ri' , 
r2',  and  r3'  and  the  corresponding  time  adjustments  to  TD3,  TD2f 
TD3,  based  on  the  geometry  presented  in  this  section. 

2.  Mathematical  Derivation  of  tlie  Oorrecticn  Factor 


An  important  question  in  this  analysis  is  hew  to  weic^t  the 
relative  precision  of  each  LDP  in  such  a  way  that  it  can  be  used  in 
the  calculation  of  a  correction  factor.  Reference  16  shows  that  for 
this  type  of  data  the  "weights  are  inversely  proportional  to  variance" 
(Ref.l6:p.60) ,  meaning  that  the  lengths  of  the  rj^'  corrections  should 
be  proportional  to  the  variances  (in  units  of  length)  of  the 
corresponding  TD  values  for  each  LDP.  Using  units  of  length,  in  the 
general  solution  for  all  sets  of  data,  the  following  equation  applies: 

r^'Li  =  [IV-6] 

where: 

rj^'  =  corrections  for  TD's  1,  2,  and  3  (in  ^isec) . 
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=  variances  of  TD's  1,  2,  and  3  (in  fjsecP-) . 
k  =  the  Oorrecticsi  Factor  (units  are  meters) 

Li  =  lanewidths  of  TD's  1,  2,  euid  3  (in  v\/iisec) , 

Note  that  <7^21^2  ig  ^  unitless  wei^ting  factor. 

For  any  triangle  there  will  be  seme  unique  value  of  k  which 
relates  the  distances  frem  point  p'  to  each  LOP  to  the  variances  of  the 
IDP's  according  to  equatic^is  IV-6.  To  obtain  a  general  egression  for 
k  vhich  can  be  easily  ccilculated  by  cenputer,  it  will  be  helpful  to  use 
only  the  knewn  values  of  the  TD's,  variances,  azimuths,  eind  lanewidths 
for  each  LOP.  The  r  values  are  readily  obtained  from  each  set  of  TD 
data  by  subtracting  the  long  term  average  vedues  from  current  TD 
Vcilues.  Vadues  of  and  are  ocxist£UTts  for  this  snail  eurea. 

If  each  LOP  is  considered  a  strai^t  line  and  a  cartesian 
oooniinate  system  is  assumed  with  north  as  the  +y  and  east  as  the  +x 
direction  of  the  coordinate  ctxes,  then  each  LOP  Cc»n  be  ^aecified  by  the 
general  equation  for  a  strai^t  line; 

y  =  mix  +  bi  (i  =  1,  2,  or  3)  [IV-7] 

Three  such  equations  can  be  written  (one  for  each  LOP) ,  giving 
three  linear  equations  in  two  unknowns,  x  and  y.  These  equations  will 
not  have  a  common  solution  vhen  the  three  IDP's  do  not  intersect  at  a 
common  point. 

For  each  master-secondary  pair,  all  IDP's  of  the  set  will  have 
the  same  slope  (mj^)  within  the  small  area  being  ccMisidered  (TDj^  is 
parallel  to  TDi^y  for  any  TD^) .  The  y-interc^3ts  (the  bj^'s)  will  be 
zero  for  the  average  IDP's  intersecting  at  point  P.  Y- intercepts  for 
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other  IDP's  with  non-zero  r  values  can  be  cedculated  from  the 


relation^ips  between  aziinuths  (denoted  a_£)  and  the  r  values  of  each 
IDP  by  the  following  equation: 


bi  =  - 

oos(aj[^-jr/2) 


sin  ctj^ 


(i  =  1,  2,  or  3)  [IV-8] 


N  (y) 


FIGURE  IV-4.  Typiced  average  LDP  through  the  origin  of  a  cartesian 
coordinate  system  cuid  corresponding  pcudllel  LDP  of  an  individucil 
TD  reading.  Slopes  are  the  scurte  for  both.  Y  interc^jt  is  the  b 
value  in  equation  for  the  IDP  (zero  for  the  average  IDP) . 
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Slopes  of  the  TD's  can  be  expressed  in  terms  of  <p  as: 


mj^  =  -tan  <p^  [IV-8b] 

Figure  IV-4  shows  the  2issunied  coordinate  system  ed.ong  with  the 
lOP's  for  which  the  slopes  and  intercepts  are  described.  When 
correcticvis  are  allied,  the  corrected  r  values  are  as  follows: 

^icLi  =  -  toi2Li2 

Note  that  in  the  equation  for  corrections  (IV-9)  the  units  of 
ricLi  are  meters.  If  is  coTsidered  a  unitless  wei^ting  factor, 

then  the  vmits  of  k  must  be  meters.  Ihis  eigrees  with  the  units  of 
EquaticMi  13. 

The  value  of  r^  is  calculated  by: 

ri  =  TDi  -  TOiav  [IV-10] 


v^iere: 

r^  =  deviation  of  an  DDP  from  the  Icxig  term  average  (Msec) . 

TD^  =  time  difference  reading  of  a  master-secondary  pair  at 
point  P. 

TDigv  =  long  term  average  of  the  TD  for  the  same  pair  at 
point  P. 

Figure  IV-3  shews  that  the  r^'  corrections  must  be  subtracted 
from  rj^  and  r3,  and  added  to  r2,  for  the  COTifiguration  shown.  It  can 
be  demonstrated  that  the  r^'  corrections  are  always  of  the  same  sign 
(both  +  or  both  -)  for  the  two  most  widely  spaced  LDP's  and  of  the 
opposite  sign  for  the  intermediate  LDP,  vhere  the  direction  vectors  at 
point  P  are  in  the  positive  directicais  of  increase  of  TO  and  the  angle 
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is  in  absolute  terms  (maxinum  being  180°).  This  is  true  for  all 
positive  and  negative  r  values.  Positive  and  negative  k  values  will 
reverse  the  signs  of  all  corrections  in  equations  IV-9  so  that  they 
will  be  correct  for  all  ocxifiguzaticns. 

The  corrected  r  values  of  equations  IV-9  cem  be  substituted 
into  equations  IV-8,  providing  a  set  of  y-intercepts  (the  b^  vzdues  in 
the  set  of  general  equations  IV-7)  for  the  corrected  IDP's  \Aiich 
intersect  at  point  P' .  Ihese  intero^±s  contain  the  unknown  Correction 
Factor  k.  Note  that  the  slopes  of  the  corrected  IDP's  are  equal  to 
those  of  the  corresponding  uncorrected  IDP's,  meaning  that  the  m^ 
values  in  equations  IV-7  cure  unchcinged. 

Substituting  the  corrected  y-interc^t  (b^)  values  into 
equations  IV-7  results  in  three  linear  equations  in  three  unknowns;  x, 
y,  and  k: 


r^Li  -  kOi^Li^ 

y  =  (-tan  02)^  + -  [IV-lla] 

CX3S  02 


r2L2  + 

y  =  (-tan  02)x  + -  [IV-llb] 

cos  02 


y  =  (-tan  03)  x  +  — -  [IV-llc] 

cos  03 
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Equations  IV-11  can  be  rewritten  as  the  follcwing  set  of  linear 


equation  in  the  unknowns  x,  y,  and  k: 


ai2i^2  nLi 

(tan  0i)x  +  y  +  ( - )k  =  - 

cos  (l>i  cos  (pi 


[IV-12a] 


(tan  <f>2)^  +  y 


(- 


-)k  = 


cos  <f>2 


^2^ 

-  [IV-12b] 

cos  <^2 


(tan  03) X  +  y  + 


( - )k  = 


cos  03 


^2^ 
cos  03 


[IV-12C] 


Solving  equations  IV-12  by  Gaussian  elimination  gives  the  follcving 
general  es^ression  for  the  Correction  Factor  k: 


Liri 


k  = 


tan  03  “•  tan  02 
cos  03 


cos  02 


b32,i2 


tan  03  -  tcm  02 


cos  03 


L2^‘^2 


tan  03  -  tan  03 

,2 -  +  L32£,^Z 

cos  02 


tan  02  -  tan 


cos  03 


[IV-13] 


vAiere: 

r^  =  differences  frcsn  average  TD's  {nsec) . 
k  =  Correction  Factor  for  three  master  secondary  pairs  (m) . 

03  =  directions  of  increasing  ID  (azimuth  from  north) . 

032  =  variances  of  IDP's  (part  of  unitless  wei^ting  factor) . 

L3  =  lanewidths  (in  meters/Msec) 

Note  that  1^32032  terms  in  denominator  are  unitless  weighting  factors. 
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Equation  IV-13  may  be  sinplified  by  the  following  substituticais: 


Ai  =  L]:  (tan  03  -  tan  03)/ (cos  03)  [IV-14a] 

A2  =  l^Ctan  03  -  tan  0i)/(oos  03)  [IV-14b] 

A3  =  L3(tan  02  -  tan  0i)/(cx3S  03)  [IV-14c] 

A4  =  +  l2CJ2^A2  +  L3CT3^A3  [IV-14d] 

Resulting  in  the  following  equaticxi  for  k: 

k  =  (A^ri  -  A2r2  +  A3r3)/A4  [IV-15] 


The  X  and  y  ooonilnates  of  point  P'  can  be  calculated  by: 

rj^cos  02  -  ^2005  0j^  ~  kuj^cos  0]^  -  ka2^cos  02 

X  =  -  [IV-16] 

(cos  0^)  (cos  02)  (tan  02  -  tan  03) 

y  =  -  (tan  a3)x  +  (r^  -  a22k)/(cos  02)  [IV-17] 

Two  other  useful  parameters  Ccin  be  ceilculated  from  the  values 
of  x  and  y.  These  are  the  azimuth  and  range  of  the  of  point  P'  fron 
point  P.  These  are  given  by  the  following : 

Azimuth  (from  north) : 

appi  =  n/2  -  tan“^(y/x)  (x  >  0)  [IV-18a] 

OpPi  =  37r/2  -  tan“^(y/x)  (x  <  0)  [IV-18b] 

Note:  vhen  x=0  the  azimuth  is  zero  for  positive  y  and  tt  for  negative  y 

Range(P  to  P'):  Rppi  =  (x^  +  [IV-19] 
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Using  the  ocnfxited  value  of  k  for  a  set  of  ooxxirrent  TO  data, 
a  correction  can  be  made  to  the  TO  vzilues  using  equations  IV-9.  Since 
IV-9  is  in  meters  it  must  be  converted  to  nsec  using  equation  3V-3b. 
The  resulting  corrected  TO's  will  intersect  at  the  point  P' . 

Equations  IV-15,  IV-16,  and  IV-17  do  meet  the  requirements 
stated  at  the  beginning  of  this  secticxi.  The  inputs  are  limited  to 
azimuths,  variancjes,  long  term  average  TO  values,  euid  concurrent  TO 
data,  and  the  values  of  k,  x,  y,  and  the  corrected  TO  values  can  be 
ocnputed  directly  from  these  by  equations  vAiich  can  be  calculated 
easily  by  a  computer. 

3.  Ncn-Dlfferential  CXsrrecticn  for  Biree  Seooniary  Stations 

When  using  the  method  of  Section  IV.  B.  2  with  a  single  Loran-C 
receiver  (i.e.  not  using  differential  methods)  a  cxtrrecticai  can  be 
calculated  based  on  existing  computer  programs.  While  not  used  in  this 
ejperiment,  the  method  of  correctioi  will  be  outlined  in  this  secti 
Since  variaticai  from  a  known  point  cannot  be  used,  the  two 
IDP's  can  be  used  as  the  basis  for  this  calculation.  Wiile  not  as 
cis  using  a  kncwn  fixed  survey  point,  this  will  provide  a  correction 
bcised  on  three  lOP's.  Ihe  method  is  outlined  belcw: 

1.  Using  the  two  LOP's  with  the  best  fix  accuracy  (available  from 

existing  computer  programs) ,  calculate  the  geographic 
position  of  the  point  of  intersection.  This  corresponds  to 
point  P. 

2.  From  the  geographic  position  of  P,  calculate  the  TO  of  the 

third  master-secx)ndcuy  pair  (the  one  not  used  in  the  previous 
st^) . 

(continued  on  next  page) 


43 


i  I  I 


3.  Using  point  P  calculate  the  azimuths  (frxxn  north)  of  the  LDP's. 

Number  the  TD's  cis  described  in  the  previous  sections,  with 
the  intermediate  TO  as  TO2  and  the  TD's  with  the  greatest 
angular  sepeuration  cis  TOj^  amd  TO3. 

4.  Use  TO  vcirianoes  calculated  over  short  time  periods  in  the 

general  eirea,  or  use  the  values  of  ineters/microsecond  cal¬ 
culated  using  equation  IV-3  as  the  af»roximate  wei^ts  of 
each  IDP.  Substitute  these  for  the  values  in  equation 
IV-13. 

5.  Assume  r  values  of  zero  for  the  LDP's  with  the  best  fix 

accuracy  used  in  step  one. 

6.  Calculate  the  difference  between  the  observed  third  TO  and  the 

calculated  value  <±>tained  in  st^  2.  Use  this  for  the  third 
r  value. 

7.  Calculate  corrected  r  values  as  in  the  previous  sections. 

8.  Using  any  two  corrected  TD's,  calculate  geographic  position. 

This  is  the  position  corrected  for  three  secondary  stations. 

Note  that  all  of  the  above  calculations  are  based  on  computer 
methods  and  programs  already  in  use  (Ref.  2)  allowing  corrected  posit¬ 
ions  to  be  calculated  readily  using  the  method  outlined. 


C.  Differential  Oorrections 


One  of  the  main  objectives  of  this  eiqjeriment  is  to  apply 
differential  corrections  cfctained  at  a  known  point  to  time  difference 
data  taken  at  unknown  points.  One  problem  encountered  was  the 
availability  of  only  one  Loran-C  monitor  receiver,  making  simultaneous 
data  collection  at  two  points  iirpossible.  Differential  Loran-C  does 
not  require  continuous  i^xiating  of  the  differential  correction  for 
general  navigation  purposes  (hourly  or  daily  corrections  cure  possible) , 
but  the  need  for  the  most  accurate  corrections  possible  for  the 
experiment  required  a  cxintinuous  i^xiate  of  the  correction. 
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1.  Mathamtlcal  Description  of  the  Differential  caarrection 


Data  frxxn  the  ejqseriment  ccaisisted  of  a  series  of  readings  at 
each  point.  Data  was  recorded  at  intervals  of  5.0  seconds  for  ituch  of 
the  experiment,  althou^  other  time  intervals  were  also  used.  Data 
consisted  of  time  versus  TD  for  \ip  to  four  possible  master-seccndary 
pairs  cind  the  signal  to  noise  ratios  for  each,  in  the  following  format: 

Time  (OW)  GRI  TDi  TD2  TD3  TD4  S/Ni  S/N2  S/N3  S/N4 

Only  three  master-secondciry  pairs  were  receiveible  in  the  area 
of  the  experment.  One  survey  point  (Range  7)  was  designated  the 
"kncwn"  point,  from  vhich  differential  corrections  were  calculated  for 
data  ta3cen  at  the  other  points.  Over  a  period  of  several  days,  data 
was  collected  at  the  kncwn  point  for  periods  of  20  to  30  minutes  at 
irregular  intervals  spaced  at  about  one  to  three  hours  (data  was 
collected  at  the  other  points  in  the  intervening  periods  as  described 
in  Section  A. 2) . 

To  provide  estimates  of  the  continuous  TD  readings  at  the  fixed 
point  during  all  time  periods  (including  the  periods  between  data 
collection  at  that  point) ,  a  least  squcires  curve  fitting  process  is 
used.  The  result  is  three  equations:  TD  as  a  function  of  time,  for 
each  master-secondary  pair. 

These  three  equations  can  be  used  to  provide  an  estimate  of  the 
TD  values  at  the  fixed  point  at  any  time.  The  difference  between  the 
calculated  TD  at  a  certain  time  eind  the  long  term  average  TD  is  the 
Differential  Correction  at  that  time.  The  ASF  corrections  at  the  fixed 
point  and  an  estimate  of  the  ASF  correction  at  the  unknown  point  (as 


cSescribed  in  Section  II-B)  also  enter  into  the  Differential  Oorrecticxi 
at  the  vmkncwn  point.  An  equaticsi  is  voritten  for  each  IJDP; 

qj  =  TOc  -  +  (ASFu  -  ASFf)  [IV-20] 

vhere: 

qj  =  Differentieil  Oorrecticxi  at  the  vmknown  point. 

TD^  =  TD  (at  known  point)  cedculated  from  regression  equation. 

~  lixig  term  average  TD  at  the  known  point. 

ASF  =  ASF  correction  at  the  known  and  lanknown  points. 

2.  least  Squcires  Regression  Equations  for  Tine  Differenoe 

Linear  regression  equaticxis  for  each  TD  can  be  written  in  the 
following  form  (Ref. 17): 


TO  =  bo  +  bit  [IV-21] 

The  bi  cxjefficients  of  the  regression  equation  can  be  found  by 
solving  the  following  set  of  linear  equations: 


nbo  +  biSti  =  ETDi 

boSti  +  biEti^  =  2tiTOi 


[IV-22a] 


An  edtemate  method  of  calculating  the  b^  values  directly  is: 


nStiTOi  -  (Eti)  (ITOi) 

bi  =  - 

nEti2  -  (Eti)2 

ZTOi  2TOi 

bo  =  -  -  bi  - 

n  n 


[IV-22b] 
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In  equations  IV-22,  etc.,  refers  to  the  sum  of 
the  time  values,  the  sum  of  the  products  of  time  and  time  difference 
values,  etc. ,  for  the  v^ole  set  of  data  available  at  the  fixed  point. 
Ihe  toted  number  of  sets  of  ID  vs  t  data  is  n. 
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V.  RESUlilS  AND  CENCIUSICNS 


Part  A  of  this  chapter  presents  the  results  of  calculations  based  on 
the  theoretical  principles  previously  discussed  and  data  obtained  from  the 
experiment.  Part  B  provides  a  summary  of  the  results  obtained,  how  these 
results  differ  from  what  would  be  expected  based  on  the  theory,  possible 
reasons  for  any  variations  from  expected  outcome,  and  a  general  conclusion. 

A.  Calculations 

1.  Azinuths,  Directiors  of  Hypei±olas,  and  Directians  of  TO's 

Azimuths  of  the  Loran-C  stations  fran  the  survey  points  are  one 
of  the  basic  requirements  of  several  calculations  used  in  this 
experiment.  Station  positions  of  Range-7  and  Luces  Point  (the  primary 
stations  used  in  this  ejqjeriment)  and  other  points  were  obtained  frcm 
Horizontal  Control  Data  published  by  the  U.S.  Dept,  of  Ccsnmerce,  Coast 
and  Geodetic  Survey  (see  ^^ipendix  A) .  Station  positions  of  the  Loran-C 
transmitters  were  obtained  frcm  loran-C  Data  Sheets  provided  by  the 
Defense  Mapping  Agency,  Hydrographic/Tppographic  Center  (see  Appendix 
B) .  For  convenience  in  avoiding  datum  conversion,  the  North  American 
Datum  1927  (NAD  27)  was  used  because  published  positions  based  on  this 
datum  were  available  for  all  of  the  station  positions  used  in  the 
experiment.  Station  positions  arc  summarized  in  Table  V-1. 
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Azunuths:  Azijmuths  fran  Range-7  and  luces  Point  to  Loran-C 
stations  9940^,  W,  X,  and  Y  were  calculated  using  a  Fortran  77 
subroutine  provided  by  CMAHTC  (Appendix  C).  Operation  of  the. 
program  was  checked  using  kncwn  data  points  eind  azimuths.  Results  of 
these  ocnputations  eire  presented  in  tabular  form  in  Table  V-2.  For 
convenience  in  visualizing  the  angular  relation-  ships,  results  eire 
also  presented  in  graphic  form  in  Figures  V-la  and  V-2a. 

Directions  of  Hyperbolas;  Ihe  directions  of  hyperbolas  for 
9940-W,  X,  and  Y  were  calculated  using  equations  IV-4  and  the  azimuths 
obtained  as  described  in  the  previous  paragraph.  A  saitple  calculaticm 
follows; 


HypertxDla  Direction  for  9940-X  at  Range-7: 


HDR7 


(Azimuth  of  9940-X)  +  (Azimuth  of  9940-M) 


2 

343.071686  +  38.358551 
2 


190.715118* 

(frcjii  North) 


Decimal  degrees  are  used  in  the  calculations  for  convenience. 
Results  of  the  calculations  cure  again  presented  in  Figures  V-lb  and  2b. 

Direction  of  Increase  of  TD;  Directions  of  increasing  TD  are 
at  ri(^t  angles  to  the  hyperbola  directions  eind  increase  in  the 
direction  of  hyperbolas  vhich  approach  the  master  station.  This 
involves  either  adding  or  subtracting  90*  from  the  direction  of  the 
hyperbolas.  For  checking  purposes,  a  Loran-C  chart  of  the  area  was 
used.  Results  cure  again  presented  in  Figures  V-1  and  2. 
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TABLE  V-1.  Statioi  Positicxis  (NAD  27) 


Station 

Latitude 

Longitude 

Range  7 

36°  39'  02.47787*'  N 

121°  49*  08.58202"  W 

Luces  Point 

36°  38’  10.524"  N 

121°  55*  38.399"  W 

9940-M 

39°  33*  07.046"  N 

118°  49*  52.241"  W 

9940-W 

47°  03*  48.594"  N 

119°  49*  52.241"  W 

9940-X 

38°  46*  57.472"  N 

122°  29*  40.050"  W 

9940-Y 

35°  19*  18.342"  N 

114°  48*  13.946"  W 

TABLE  V-2.  Azimuths  Fran  Rcmge-7  and  Luces  Point. 


Range-7  to 

Degrees ,  min ,  sec . 

Decimal  Degrees 

9940-M 

38°  21*  30.784" 

38.358551° 

9940-W 

7°  32*  22.823" 

7.539673° 

9940-X 

343°  04*  18.070" 

343.071686° 

9940-Y 

101°  47*  36.229" 

101.793397° 

Luces  Point  to 

9940^ 

38°  51*  08.262" 

38.852295° 

9940-W 

7°  45*  43.837" 

7.762177° 

9940-X 

344°  23*  49.981" 

344.397217° 

9940-Y 

101°  36*  16.319" 

101.604533° 
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9940-W  7.539673 


FIGURE 


9940-Y 


FIGURE 


FIGURE 


9940-X  343.071686 


V-la.  Azinuths  Fran  Range-7  to  9940-M,W.X,and  Y. 


^-Ib.  Directiais  of  Hypertxslas  at  Range-7. 


9940-Y 


9940-W  112. 949112* 

^-Ic.  Dinecticns  of  Maxinum  Increase  (Gradient)  of  ID's  at  R-7. 
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FIGURE 


9940-Y 


FIGURE 


FIGURE 


9940-X  344.397217 


V-2a.  Azimrths  Fran  Tuoes  Paint  to  9940-fl,W.X,cUTd  Y. 


V-2b.  Directions  of  Hyperbolas  at  luces  Point. 


9940-Y 


9940-W  113.307236* 

^-2c.  Directions  of  Maxinum  Increcise  (Gradient)  of  TO's  at  IP. 
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2.  lanewidths 


Lanewidths  were  calculated  frtan  Equation  IV-3a  using  the 
azimuths  oonputed  in  Section  V-A-1.  Note  that  the  units  of  lanewidth 
ctre  ra/^lsec.  With  a  standard  leme  of  1  fisec  the  width  of  a  lane  is  the 
same  number  in  meters.  Velocity  of  transmission  was  calculated  using  c 
emd  the  index  of  refraction  specified  in  the  table  of  loran-C  constcints 

supplied  by  the  Defense  Mapping  Agency  (Ref.  4).  Sample  calculations  of 
lanewidth  follow: 


c 

Velocity  of  Transmission  =  -  [V-l] 

n 

2.99792453  x  lO^  m/sec 


1.000338 

=  2.996911624  X  10®  Vsec 

lanewidth  of  9940-Y  at  Range-7  for  standard  1  ijsec  lane: 
(10"®sec)  (2.996911624  X  10®  itv/sec) 

1^7  = 

101.793397*  -  38.358551* 


=  285.024  m 


Results  of  the  calculation  of  lanewidths  are  given  as  the  L 
values  in  Table  V-3.  In  the  table,  L^,  L^,  and  L3  correspond  to  the 
Icinewidths  for  9940-W,  X,  emd  Y,  respectively. 
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3.  Asrerage  TD's,  Standard  Deviations,  Varianoes 

Because  of  the  leirge  amount  of  data  involved,  it  was  neoessciry 
to  compute  statistical  information  losing  a  computer  program.  Since 
data  was  in  the  form  of  three  TO's  (one  for  each  master-secondary  pair) 
every  five  seconds,  and  readings  were  taken  for  periods  of  20  to  40 
minutes  at  each  station,  a  Fortran  77  program  with  subroutines  for 
statistical  information  weis  written  (see  i^Y>®rdix  C) . 

Station  Pange-7  was  the  ''known''  station  of  the  differential 
system,  at  which  more  readings  were  taken  than  any  other.  For  this 
reason,  and  because  regression  equations  of  TD  vs  Time  were  needed  only 
at  that  station,  statistical  information  was  calculated  only  for 
station  Range-7. 

Mean  of  the  TD's  was  computed  by  the  following  (Ref.  18): 

CTDi 

Mean  TD  =  -  [V-2] 

n 

Variances  of  the  TD's  were  calculated  using  the  following  (Ref.  18): 

_  (j7rD|)2/^ 

s2  =  jV-3] 

N  -  1 

Standard  deviations  were  calculated  from  variances  by  (Ref.  18): 
formula  (Ref. 18): 

s  =  (s^)^  [v-4] 
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Equations  V-2,  V-3,  and  V-4  v^ere  incorporated  into  a  Fortran  77 
ocnputer  program  vAiich  was  used  to  calculate  statistics  of  TD's  for 
9940-W,  X,  and  Y. 

It  should  be  noted  that  these  statistics  were  in  terms  of  time 
differences,  \diile  in  seme  ca5=yi.s  it  was  necessary  to  know  standard 
deviations  and  variances  in  terms  of  distances.  Equations  IV-5  were 
used  for  this  conversion  vAiere  necesscury. 

Corrected  vs  Uncorrected  Data;  One  use  of  the  statistics  of 
the  data  was  to  cenpare  the  stcindard  deviation  of  raw  data  to  the 
standard  deviation  of  data  vrtiich  was  corrected  for  three  secondary 
stations.  It  Wcis  considered  unlikely  that  the  three  station  correction 
could  remove  bicis  from  the  TD  readings  (error  in  terms  of  absolute 
position)  because  most  of  this  is  due  to  systematic  error  rather  than 
randan  variations.  It  Wcis  believed  that  the  three  station  correction 
would  be  most  effective  in  reducing  the  effect  of  randan  variations  in 
time  differences.  The  statistics  on  uncorrected  and  corrected  data  show 
results  consistent  with  this  expectation  (Table  V-7). 


4.  Error  Ellipse  and  Circular  Error  at  Range-7 

Range-7  was  used  eis  the  known  point  of  the  differential  system, 
and  a  larger  amount  of  data  was  collected  at  that  point.  Using  the 
crossing  angles,  standard  deviations  and  variances  of  9940-X  and  Y  (the 
best  pair  at  that  point) ,  dimensions  and  orientation  of  the  standcird 
error  ellipse  and  CEP  50%  and  90%  circle  radii  were  calculated  for  the 
point  using  the  information  outlined  in  section  III-B. 
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The  standard  error  ellipse  and  emor  circles  for  50%  and  90% 
were  plotted  over  a  scatter  plot  of  data  from  Range-7,  ihe  result  is 
shown  in  Figure  V-3. 

Calculaticxis  of  the  standard  error  ellipse  dimensions  were  as 

follows: 


Using  9940-X,Y: 


Cl  =  O9940_y  =  (0.072546) (285.024)  =  20.677m 

C2  ~  ^9940— X  ~  (0.068638) (322.024)  =  22.103  m 

0  =  250.075974”  -  190.715118°  =  59.360856” 


Calculating  the  semi-major  and  semi-minor  axes  of  the  ellipse: 


2  _ 


2  _ 


20 . 6772+22 . 1032+  ( (20 . 677^+22 . 1032)  2-4  (20 . 677^^2 .  lO32sin20) ) 

2sin20 

30.595  m 

20.6772+22 . 1032- ( (20. 6772+22.1032) 2-4 (20. 677^^2 . lO32sin20) ) ^ 


2sir)20 


=  17.362  m 


Orientation  of  the  major  eods  of  the  ellipse: 


tan  20  = 


20.6772sin  2(59.360856”) 


20.6772  cxs  2(59.360856”)  +  22.1032 
0  =  26.474268” 


Azimuth  of  major  axis  =  250.075974”  -  26.474268  =  223.601716” 

=  223.601716”  -  180”  =  43.601716”  (frcsn  North) 
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Using  tables  in  Reference  6,  pp.  18-20,  radii  of  the  50%  and 
90%  error  circles  were  calculated  ais: 

CEPgQ  =  27 >936  in  CEP^q  =  63.226  in 

Because  the  data  inclucSed  in  the  plot  of  Figure  V-3  includes 
two  intervals  of  data  collectiai,  the  points  are  not  as  evenly 
distributed  as  they  would  otherwise  be.  For  this  reason  there  are 
fewer  points  near  the  center  of  the  plot  than  would  normally  be  the 
case.  The  error  circles  eind  ellipse  seem  to  be  reasonably  correct  for 
the  given  data. 

5.  Garrection  Factxir  for  Three  Secondary  Stations 

The  correction  for  three  secondary  stations  was  designed  as  a 
technique  to  reduce  the  randcm  error  associated  with  time  difference 
data.  It  was  e}g)ected  to  increase  the  precision  of  the  position  rather 
than  the  absolute  position. 

For  a  variety  of  reascns  there  will  be  sane  covariance  between 
concurrent  time  difference  readings  from  different  master-seocaidary 
pairs.  The  correction  for  three  secondcury  stations  is  not  expected  to 
reduce  the  error  associated  with  this  type  of  random  error.  The  nai- 
oorrelated  random  error  in  ID  readings  should  be  reduced  by  this 
correction  to  some  extent. 

In  order  to  reduce  the  nuntoer  of  cnpitations  and  the  time 
required  to  process  data,  most  cOTstant  factors  eissociated  with  the 
correction  factor  can  be  corputed  in  advance  and  used  as  constants  in 
the  conputer  program.  This  reduces  the  correction  to  only  a  few  lines 
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of  relatively  siirple  examputation.  Because  of  the  letrge  amount  of  data, 
this  reduces  significantly  the  time  required  to  cerpute  corrections. 

In  the  calculations  of  adjustments  to  positions  of  Range-7  and 
Luces  Point,  the  only  factors  which  differ  significantly  cure  the  long 
term  time  differences  at  the  two  points.  In  this  case  the  same 
factors  can  be  used  for  correcting  both  positions  for  three  secondary 
stations.  In  cases  where  seme  of  the  parameters  may  differ 
significantly,  it  would  be  necesscury  to  calculate  different  factors 
for  each  position.  Table  V-3  gives  parameters  for  irput  to  the  Fortran 
77  correction  program  for  calculation  of  the  Aj^  parameters  and 
calculation  of  the  corrections  to  data. 

TABLE  V-3.  Parameters  used  in  calculation  of  correction 
factor  k . 


Parameter 

Value 

<^1 

112.949112' 

100.715118' 

<p2 

340.075974' 

0.150606  Msec 

^2 

0.068638  Msec 

03 

0.072546  Msec 

Ll 

563.934  m/nsec 

L2 

322.966  nv/Msec 

285.024  nv/Msec 
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6.  Regression  Bgiiatlcns  for  Time  Difference  vs  Time 


Since  only  one  lorcin-C  monitor  receiver  was  available,  it  was 
necessary  bo  develop  a  method  of  estimating  TD's  at  the  known  point 
during  times  vtien  differential  corrections  were  required  at  the  unknown 
points.  Linecir  iregression  equations  of  TD  vs  time  were  used  for  this 
purpose. 

Using  equations  IV-21  eind  IV-22,  the  data  from  Range-7  Wcis 
processed  to  obtain  the  required  sunmations  of  time  and  time  difference 
necessary  to  calculate  the  bo  and  b]^  coefficients  of  the  linear 
regression  equation.  Only  Range-7  data  from  the  nearest  time  periods 
before  and  after  the  differential  correction  period  for  Luces  Point 
data  were  used. 

The  procedure  for  processing  the  data  was  first  to  apply  the 
correction  for  three  secondary  stations  to  the  Range-7  and  Luces  Point 
data  during  the  time  periods  of  interest.  Range-7  data  was  then 
processed  to  obtain  the  necesseiry  summations  for  the  calculation  of  the 
b  parameters  of  the  regression  equations.  Using  the  regression 
equations  for  Range-7  data  and  long  term  TD  averages  at  the  two  points, 
differential  corrections  for  the  Luces  Point  data  were  calculated  as 
described  (Section  V-7).  Fortran  subroutines  (Appendix  C)  were  added 
added  to  the  statistics  program  to  obtain  the  summations  and  a  separate 
program  was  used  for  the  differential  correction  of  luces  Point  data. 
Results  of  the  corrections  were  ccsipared  to  calculated  values  of  TD's 
provided  by  the  Defense  Mapping  Agency,  Hydrographic/Topographic 
Center.  TD  vs  time  is  plotted  in  Figures  V-4  along  with  the  regression 
line  for  the  data.  Regression  equations  are  in  Table  V-4. 
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TD  (uaoc) 


16316.9. 


TIME  vs  TD 

9g40-W  9/28/87 


16316.4. 


16316.1 


TIME  (PST) 


Figure  V-4a .  TD  vs  time  and  the  regression  line  for  9940-W 
for  Range-7  data  corrected  for  three  secondary  stations. 
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S.O  iO.O  il.D  i2.Q  i3.0  14.0 

TIME  (PST) 

Figure  V-4b.  TD  vs  time  and  the  regression  line  for  9940-X 
for  Range-7  data  corrected  for  three  secondary  stations. 
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Figure  V-4c.  TD  vs  time  and  the  regression  line  for  9940-Y 
for  Range-7  data  corrected  for  three  secondary  stations. 


TABTE  V-4.  Regression  Equations  of  TD  vs  Tiitie  at  Rcinge-7. 


9940-W: 

TD^  = 

16316.51564 

-  0.00807243t 

994 0-X: 

TDx  = 

27522.33725 

-  0.01385199t 

994 O-Y; 

TDy  = 

42749.29371 

+  0. 010567 07t 

time  units:  hours  with  minutes  and  seconds  in  decimal  hours. 

7 .  Differential  Corrections 

Equation  IV-20  describes  the  differential  correction  vrtiich  is 
applied  to  TD  readings  at  an  urOcncwn  point  based  on  TD's  at  a  known 
point.  In  this  experinvent  the  positions  of  all  points  are  known,  based 
on  their  published  positions  from  survey  data  (see  appendices) .  One 
point,  Range-7,  was  chosen  eis  the  known  pjoint  for  pxjrpxjses  of 
differential  corrections.  Data  from  one  other  point  was  chosen  as  the 
unknown  point  for  purposes  of  ccsiparison  to  its  known  position. 

Ihree  things  must  be  known  to  calculate  a  differential 
correction  based  on  equation  IV-20.  Rather  than  use  long  term  averages 
in  this  equation,  the  expected  TD's  at  both  points  were  corputed  by  a 
standard  Loran-C  program.  Current  TD's  at  Range-7  were  calculated  from 
the  regression  equations  described  in  the  last  section.  The  ASF 
corrections  in  equation  IV-20  were  calculated  for  each  point  by  the 
Defense  Mapping  Agency  using  a  computer  program  which  is  used  in  the 
preparation  of  ASF  correction  tables,  but  which  can  be  used  for 
determining  the  correction  at  individual  points.  The  difference  in  ASF 


64 


corrections  is  used  in  equatic^i  IV“20.  Table  V~5  gives  expected  TO 
values  (used  as  Icxig  term  averages)  at  the  two  points  of  interest: 

TARTK  v-5.  Calculated  TO's  used  as  Icxig  term  averages  (in  ^sec) . 


Range-7 

Luces  Point 

9940-W 

16314.68 

16300.09 

9940-X 

27523.13 

27496.79 

9940-Y 

42749.49 

42755.87 

ASF  corrections  were  approximately  the  same  at  both  points  due 
to  the  nearness  of  the  points.  Differences  in  ASF  corrections  used  in 
equation  IV-20  were  very  snail:  0.005,  0.002,  and  0.003  /xsec  for  9940- 
W,  X,  and  Y,  respectively.  Corrections  of  this  magnitude  are  probably 
less  than  the  acxuracy  of  the  data  used  to  calculate  them,  but  at 
longer  ranges  these  may  be  large  enoui^  to  be  significant. 

Ihe  final  equation  for  the  differential  corrections  based  on 
TD's  at  Range-7  eure  in  Table  V-6.  Note  that  t  is  Pacific  Standard  Time 
in  hours,  with  minutes  and  seconds  expressed  as  a  decimal. 

TABIZ  V-6.  Differential  Corrections  (C^  in  /isec,  t  in  hours  PST) 


^^idw  ~ 

-1.83064  + 

0.00807243t 

0.79575  + 

0.01385199t 

C!dy  = 

0.19829  - 

0.01056707t 
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8.  Accuracy  of  Position  Fix 


Figure  V-5a  shows  a  plot  of  unoorrected  Range-7  data  beised  on 
ej^jected  ID's  calculated  froii  the  )<nown  geographic  position.  Ihe  grid 
is  in  meters,  oriented  north-south  emd  east-west,  euid  with  the  center 
at  the  published  geographic  position  of  Range-7.  When  the  Correction 
For  Three  Secondary  Stations  is  applied  and  the  average  of  the  Range-7 
data  is  placed  at  the  origin,  the  result  is  as  shewn  in  Figure  V-5b. 

Note  that  the  plot  is  new  centered  on  the  origin  due  to  the 
fact  that  all  absolute  error  (bicis)  has  been  removed  fron  the  data. 
The  three  station  correction  has  also  decreased  the  standaird  deviation 
and  varicince  of  the  Range-7  data  as  shewn  in  Table  V-7. 

Figure  V-6a  shews  Luces  Point  data  before  application  of  any 
corrections.  These  pxiints  cure  plotted  based  on  the  expjected  TD's  at 
the  knewn  coordinates  of  the  px>int.  The  center  of  the  uncorrected 
Luces  Point  data  is  at  x  =  375.820  m,  y  =  85.501  m,  and  at  a  range  of 
385.423  m  cind  at  an  azimuth  of  77.183°  (frem  North)  frean  the  origin. 
Radii  of  the  CEP50  and  CEPgo  error  circles  are  14.9  m  and  28.4  m, 
respjectively. 

Figure  V-6b  shews  Luces  Point  data  after  the  application  of 
both  three  station  and  differential  corrections.  Again,  the  origin  is 
at  the  knewn  pxssition  of  the  pxjint.  For  visual  conparison  pxirpxDses, 
the  same  scale  is  used  for  both  Figures  V-6a  and  b.  For  the  corrected 
data,  the  center  of  the  data  is  at  x  =  47.213  m  and  y  =  -7.015  m,  with 
a  range  of  47.731  m  and  at  an  azimuth  of  98.451°  (from  North)  from  the 
origin.  CEP5Q  and  CEPgQ  radii  were  12.6  m  and  26.0  m,  respectively. 
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UNCORRECTED  9940-X,  Y 
SCATTER  DIAGRAM 
STATION  RANGE-7 
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Figure  V-5a.  Scatter  Diagram  For  Rcinge-7  Uncorrectad  Data.  Scale  is 
in  meters  oriented  N-S  cind  E-W,  centered  on  the  time 
differences  calculated  from  the  published  coordinates 
of  the  point. 
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CORRECTED  9940-X.  Y 
SCATTER  DIAGRAM 
STATION  RANGE~7 
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Figure  V-5b.  Scatter  Diagram  for  Range-7  Data  Cjorrected  for  Ihree 
Secondary  Stations  and  Differential.  Scale  is  in 
meters  oriented  N-S  and  E-W.  Center  of  data  is 
corrected  to  correspond  with  calculated  time 
differences  for  the  published  coordinates  of  the  point. 
See  text  for  details. 
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UNCORRECTED  9940-X.  Y 
SCATTER  DIAGRAM 
STATION  LUCES  POINT 


Figure  V-6a.  Scatter  Diagram  For  Tiioes  Point  Uhoorrected  Data.  Scale 
is  meters  oriented  N-S  and  E-W,  centered  on  the  time 
differences  calculated  from  the  published  coordinates 
of  the  point. 
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SCALE  IN  lETERS 


Figure  V-€b. 


Scatter  Diagram  for  luoes  Baint  Data  With  Differential 
Oorrecticns  Based  on  Range-7  Data  and  With  Oarrection 
For  Dvree  Secxjndciry  Stations.  Origin  of  plot  is  at  the 
time  differences  calculated  frcn  the  publi^ied 
coordinates  of  the  position. 
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TABLE  V-7.  St3ndaid  Deviaticn  and  varianoe  of  Range-7  Dafa  With  and 
Without  the  Ocsrrecticn  For  Three  Secondary  Stations. 


Pair 

Without  Oarrection 

With  Oarrection 

a 

0.068638 

0.051161 

9940-X 

0.004711 

0.002617 

a 

0.072546 

0.052081 

9940-Y 

0.005263 

0.002712 

B.  CSonclusicns 

Based  on  the  results  given  in  part  A,  the  absolute  accuracy  of 
the  Luces  Point  Data  was  inproved  from  eui  average  error  of  385.423 
iteters  to  an  average  error  of  47.731  meters.  Practically  all  of  this 
iitprovement  was  a  result  of  application  of  the  differential  correction. 

Precision  of  the  luces  Point  data  was  inproved  from  14.9  m  and 
28.4  m  for  the  50%  and  90%  CEP  radii,  uncorrected,  to  12.6  m  and  26.0  m 
for  the  corresponding  radii,  corrected.  Most  of  this  inprovement  is  a 
result  of  the  correction  for  three  secondciry  stations. 

These  figures  are  bcised  on  the  processing  of  486  sets  of  TD 
data  frcsm  Range-7,  the  "known"  point  of  the  differential  system,  and 
252  sets  of  TD  data  at  Luces  Point.  Each  set  of  data  included  time  cind 
three  TD's,  one  for  each  master-secondary  pair  (9940-W,  X,  and  Y) . 

Under  ideal  conditions  a  better  absolute  correction  would  be 
expected  for  these  points.  From  em  examination  of  the  characteristics 
of  the  data,  the  weakest  point  is  prchably  the  lack  of  a  continuous  set 


of  TO  readings  at  Rzmge-?.  Error  in  the  estimate  of  TO's  based  on  the 
linear  regression  lines  may  be  as  much,  or  more  tham,  the  48-meter 
average  error  in  the  cxjrrected  positions  calculated  for  Luces  Point. 

Other  sources  of  error  may  include  the  fact  that  sene  of  the 
Loran-C  sigrals  traversed  some  seawater  path  near  Luces  Point,  with 
unknewn  results  c»i  TO  readings  at  that  point.  This  type  of  error  is 
known  to  result  uncier  similar  cxsiditicxis,  but  no  means  of  estimating 
its  magnitude  is  available  at  this  time. 

In  sumnary,  the  absolute  positions  calculated  from  differential 
cx)rrec:tions  of  Loran-C  data  seem  to  be  a  significant  inprovement  cjver 
the  positions  calculated  from  raw  cLata.  While  the  carrecotion  for  three 
secondary  stations  does  not  add  much  to  the  estimate  of  positions  in 
this  experiment,  the  inprovement  of  50%  CEP  radii  by  about  15%  does 
seem  significant  encxi^  to  warrant  further  cxaisideration. 
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APPENDIX  C-1 


C  THIS  PROtRAH  READS  DATA  PROH  A  PILE  OP  TIRE.  A.  I,  C 

C  CONTINUES  TO  THE  ENO  OP  THE  PUE.  CALCULATING  THE  HEAN, 

C  VARIANCE,  AND  STANOARO  DEVIATION  OP  THE  A,  B.  AND  C  COLUHNS, 

C  AND  PRINTS  OUT  THE  RESULTS  INTO  A  PILE  tdout.f 

deublt  drtclilon  t lat . t ia*! .1 i**}, t iatS. •.b, c.nk, ni.nb, nc , 

* tua, tuat, (uabi tuac, auaaj, auabl, tuac2. auaal . auabi ,iuacl,Bl,bl,cl 
earn  (unital.fllt  a  'outt.f*.  fora  a  ‘foraattad*. 
aaccriia'aaourntlar  ) 

C  eorn  (unlta3,fll«  a  'outA.f',  fora  ■  ’foraattad’, 

C  aaccfi A* ' aaoua nt ill ‘ I 

rtwind  I 

C  raadll,*)  fritln 

C  rtadd,*)  itcln 

so  uritt  (*,Al  tia(l.tlat2,tlaaS.il,bl.cl,na.na,nb,nc 
tlBC  a  tiifl  «  tia(2/i0.  *  tiar3/S600. 

C  vrlta  (2, SI  I iat, al ,bt , el 

lul  a  AUa  *  1. 
tuaa  a  luat  «  a 
auaal  •  auaal  <  il 

auai2  a  auaal  «  a'l 

auab  a  auib  «  b 
auabi  a  auibl  «  bt 

auabi  a  auabi  «  b*b 

auar  >  auae  *  c 
auad  a  auacl  •  el 

auad  a  auacl  <  e‘e 

raadll.I.arrilO.ioatitaiofrr)  tia«l,tia«I,tiae!,al.a,bl,b, 

<c  1 ,  c 

if  (iofrrl  10. SO. 10 

10  aua  a  sut  •  l. 
ur i I « ( • .  1 1  lof rr 

eall  aliauaal. tuaa, aua, auaal  I 
call  all auabi, auab, aua, auabi 1 
call  til auacl. auae. aua, auacl I 
2  foraat Itl .dl.O. li.d:.0,t7,d2.0,tll.dS.S, tll.dS. ). tIT, 

«dS. ],t 30,dS.3, tSA.dS. S.ISV.dS. SI 
a  iroraitlfS.O,2i,fS.O,2i,fS.O,3i,flO.S,2i.flO.S,Ii,flO.S,2i, 

<lt,it.it,ial 

A  f oral! ( ' i o»rr  •  ',131 

8  roraitlfl3.8.2i,fl0.3,3i,fl0.3,3i,fl0.3l 

9  foraat ( ’data  arror' ) 

AO  atop 

«nd 

C 

tubreut i nc  al I il, i <  d, a  1 1 

doublf  oracltlon  il,  i.d.atin.varlan, ad. at 

aian  •  al/d 

varlan  a  I  iI*  ■  * i /d  I /( d- 1 1 
ad  •  daqrt(virlan) 
wriltl*,ll)  Btan, varlan, ad 
C  vrllflltll)  atm, varlan, ad 

raturn 

11  feraatCiaan  a  ' ,  f  1 2. 3, 2i .  ’  vir  1  anea  a  ' ,  f  1 1. 1,  Ii , '  ad  •  ',fl2.|} 

and 


78 


APPENDIX  C-2 


10  r«»  - - - - ...... 

:0  rta  •••  fROCRAn  FOR  (AUSSUM  ILItllNATION  UITH  RIVOTAL  CONOCNSATION 


30  rta .  FOR  SOLUTION  OF  SYSTCnS  OF  LINtAR  EOUATIONS 

40  r*» .  UITH  UF  TO  10  ONANOHNS  . 


50  re» . - . 

to  rf* 

70  torn  *1'.  II.  'Yttst.f’  :  Ola  AKlO.lDidla  BKlDidia  RI(M):dla  illll) 
to  Input  il,  rows!,  cols! 

to  r  >  rows!  :  e  •  eels!  :  rowl  ■  1  :  cell  •  1  :  hloHrow  >  1 
95  rta 

100  r»a  -- — -  Input  tha  aatrii: 

105  rca 

1 1 0  Tor  J  •  1  t 0  rows ! 

120  Tor  k  •  1  to  cols! 

130  Input  II,  AI(J,k) 

140  IT  k  <)  cols!  than  liO 
150  AKJ.kl  a  -1  ‘  Al(J.k) 

ItO  neit  k 
170  nait  J 
175  raa 

110  raa - baein  aain  proiraa - - - - - - - .............. 

115  rta 

190  tosub  500  :  rta  ---  chack  IT  unloua  solution 

200  oosub  ttO  :  raa  •••  prints  out  arrar 

210  eosub  ISO  :  rta  cnacks  Tor  row  with  htyhast  Tirst  non>iaro  antra 

21 5  raa 

220  if  highrow  •  rowl  than  240 
225  rat 

23C  ootub  700  i  ras  •••  tichanga  rows  iT  nactsssra 

240  goiut  760  !  raa  •••  sals  Tirst  non  laro  antra  in  rowl  to  1.0 

250  gosut  820  !  raa  ---  sats  Tirst  antra  to  tare  In  following  rows 

255  re» 

2L0  rowl  «  rowl  4  1  ;  coll  •  coll  ♦  1  :  highrow  a  rowl  :  raa  ..*  incrtaant  tounla'i 
245  raa 

270  if  rowl  ()  rows  than  200  :  raa  •••  checks  IT  last  row  has  bear,  rtachad 

280  gokub  760  :  raa  sats  first  non  laro  antra  in  linsl  row  to  1.0 

290  gokub  t90  :  rat  •••  final  printout  of  arrta 

295  raa 

300  raa.  -.•••  calculatas  the  unknown  i  vaiuas: 

305  rca 

310  ilicols!)  >  1  :  a  •  rows ! 

320  Tor  I  a  a  to  1  Step  -1 

330  for  1  a  1  to  cols!  •  1 

340  tSfl)  a  il(i)  .  iKlal)  •  Aid, jail 

350  nait  j 

360  nait  1 

345  raa 

370  raa  prints  out  the  i  values; 

375  rta 

380  Tor  k  a  I  to  rows! 

390  print  ’B'  k-i  '  a  •  il|kl 
400  nait  k 


410  and 

420  rca  .  and  of  aain  prpgraa 

425  rca 

430  rat  begin  subroutinas  • 

435  raa 


500  for  1  a  1  to  rows! 

510  for  3  a  i  to  rows! 

520  if  I  a  j  than  610 
530  Tor  k  a  1  to  cell! 

540  if  Aili.kl  a  0  than  160 
550  Nllkla  AIIJ.kl/Alli.k) 
560  nut  k 
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ft70  f»r  ■  ■  1  to  eel*l  •  1 
sio  if  ni(«)  ()  ni(«*i)  thtn  ilo 
StO  ntit  ■ 

»00  print  'CQUATIONS  NOT  INDEfENOCNT  ...  HO  UNIOUE  SOLUTION  ...  ONE  SOLUTION  lOLLOUS’ 

tiO  ntil  i 

620  not  1 

6S0  r*turn 

63S  rtP 

6(0  rts  .  lubroutint  to  ehtek  for  raw  with  hlphtpt  firtt  npn>(tro  tntrrs 

650  (or  i  •  rowl  to  rowd  •  I 

660  if  pbilAKl.eell))  I  ■b((AI(l4l,eoll))  thtn  610 
670  highroH  «  1  «  I 
650  ntit  1 
690  return 

695  rti  -  tunroutint  to  tiehongt  row*  pno  put  hightit  firtt  pntrr  firtt: 

700  for  1  >  cell  to  colt! 

710  BKll  •  Alirowl,  i) 

720  Alirowl, 1)  •  AKhighrow,  1) 

750  Allhighrow,  11  >  SKI) 

7(0  ntit  i 
750  rtturn 

755  rtP  •••••  lubroutint  to  itt  firtt  tntri  in  firtt  row  to  1.0: 

760  SI  >  Alirowl, coll) 

770  for  k  •  coll  to  coll! 

750  if  Alirowl, colli  •  0  thtn  510 

790  Alirowl. k)  •  Alirowl,  k)/SI 

791  print  k 
500  ntit  k 
510  rtturn 

515  rto  lubroutint  to  itt  othtr  firtt  tntritt  to  ttro: 

52C  (or  J  •  rowl  ♦  I  to  rowt! 

550  PI  •  All j. (01 1 l/Allrowl,eol 1 1 
5(C  for  k  •  coll  to  coll! 

550  AIIJ.O  >  AlU,k)  •  Alirowl, kl’PI 
560  neit  k 
570  neit  J 
550  rtturn 

555  rti  >••••  lubroutint  to  print  out  the  oitrli  A; 

590  for  J  •  1  to  rowt! 

900  I  or  k  •  1  to  colt! 

910  print  AIU,k), 

920  ntit  k 
930  print 
9(0  ntit  J 

950  print  rowl:  - ............ - ......................................... 

960  rtturn 

970  print  'NOT  ENOUGH  EQUATIONS  ...  HANY  SOLUTIONS  AAE  POSSIBLE’ 

950  rtturn 
990  tnd 
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FH06l'*n  *7" 

*  ■  (]7e;c».i 

e  • 

eFtN(uNM»i.r  lit  > '  •»  if>.  f '  ,ro»n*'re»n*T7ef‘ . 

•«CCtSS>‘tC«U(N1|Al' ) 

KtUiND  1 

C  Olout)''  It  fro*  DOint  I  to  point  7 

ltt*Ml,9l 

ltt*D(1.7l  *L*T|,AL0N| 
flAC  ■  0.0 

CALL  INVlA,B,AlA1|,Al0Nl.ALA1},AL0N}.rLAT,Dt!T.A7l 
A7  •  A7  •  3.U1S«?(S« 
a:  •  A7/«.  WlSOJASCltC.O 
ir  IA7-J60.0)  10.10.20 
20  A7  •  A7  •  JOO.O 
10  wAMtC.J)  DISt.*2 

>  •  alat;/.V  UlSO2tiA‘l(0.0 

Y  •  Aicc/J.uisorosi'uo.o 

V  •  AlATl/3.  UlS»:<:<MtC.O 
U  I  AlCM  /}.  U1  S«2e  M  *  1  tC.  0 

URntC.il  K.Y 

UFIUC.S'  V.L' 

?  rofn»T (r«, I. Ti :. r«. 1 1 

j  f ofr.iT I ‘Di jTANc t  •  ■ , n :. j. I « . *nt TtFi" . j>. 'AnnyiH  •  '.nr.e.:*. 

. • I NOFTh ' •  I 

(  rprnttCrtlA”.  •,r«/lAT:’,2».rt.3.‘N'.S«.‘L0H;'.2».f».i.‘w‘' 

ts': 

c 

toe*  COT  1  At  lAYItOf  At  .  FORAI'.HTI  .AlAl  .»ti;.»lN:.CflO.MtT.  a:1 

1'  icfic  ict.ici.io: 

ic:  tlit  .  1 . c.fC'AAt -tiRAr  • 

ta.at;  I  rt.»'’rLAi 

f  1  ■  RAI at:* 1 . 2! 

r:  <  ta.at:*:, I 

rj  .  fAtAT2*c. 
fA  .  »alat;<c..  itf 
fi  t  f ALAT2*O.Ot;! 
rt  t  fii.Ai;4rLAT 

r:  ■  r 0*1.0 
re  •  ro’o.i 
PI  I  J.lAI^02t54 

Tupri  1  o.oejiesjc* 
cric  •  1.0 

lO:  ftlAI  •  ATANI 1 1 . O-riAt I ■ JINIPlIl l/COSIAlTl 1 1 
eecTAi  •  siNietiA] ) 
ctCTA]  •  cctietTAii 

ttlAt  ■  ATAN((I.O-riAl)‘5IN(Pn2)/C0S(AlT:T) 

Sf£TA2  .  SIMetTA;) 
cetTAj  •  coJiftTA;) 

Otu  •  PlNI-PLA': 

AOCLl  >  AtilOtlU 
ir  lAOtll-PIl  IOt.lC!,10I 

ICT  Attu  •  tUOrj.AtUl 

lot  liotl  •  ilAlAOtUl 
CM'tl  •  (CSIAOtU  I 
A  •  JttTAl’fUTA: 
e  ■  ce'.TAi'cttTA; 

COtHl  •  A*f‘COtitl 
Tirnl  • 

ltOrTHilOtl'(et1A2l”:*(J(ItA2'CfrTA|.J(HAI'Clt»A2*COOtL)**.*t 
C  •  e'$10tL/»lAMl 
tfl  •  1,0'ft 
P*Tl  •  ASIMtIPMIl 


i 


81 


appendix  C-3 


ir  (corNi)  20M0i,}0i 

J07  ^MI  ■  ri-fMI 
201  PNISO  •  PMIVHI 
CJ^Ml  • 

CTPHI  •  COPHI/SIPHI 

psreo  ■  8iPMi*eoPHi 

Tinni  •  r7*PMi 

u»n2  •  *'(r4*iiPMi-f2‘PMi$o*c»PHn 

Tinns  •  tJ(Mf2*fMI10*CTfHl-fl*(PKt4flYCO)) 

T(Kn<  •  A‘*'r2*PSTC0 

itiins  «  ifi'thMr5*(PHi4P$Yeol-r2*fHiJO*ctPHi-r4*PiTCO*coPHi*cofHi  i 
YCPn6  ■  **IB*P2*(PHl$0*ClPHl«P»YeO*COPHl) 

D!$T  •  P0P*0*(TCPni«1(Rn24UPRS-1tPn4*Te*ns«Tritni) 

YERn7  •  P4*PHI 

TC»m  •  **lr2*$lPHl4f*tAT2*PMI$0*C»PHt) 

tt»B4  ■  IBMf J*r J*P$YC04r*t.*Y2*PHI$0*CtPHI-rj*PHn 

2LAn  •  C*(TCRn7-TCAnt4T{Ant)4A0tt.L 

CTA24($(CTA2‘CaETAl-C0S(2LAn2)‘SttIAl*CttTA2)/(SIN(2LAn)*ceCTA2) 

IP  (CtA2)  210,204,210 
204  C7A2  •  O.OOOOOOOS 

210  A2  •  ATAN(l.0/eiA2) 

IP  (OCLL)  215,214, 211 

211  IP  (DPLL'PII  211.212,212 

211  IP  IC1A2I  220,221.221 

220  A2  •  A24P1 
SOlO  221 

212  IP  (CTA2)  217,216,212 

215  IP  (DCLliPlI  211,211.214 
214  IP  (C7A2I  217,214.216 
217  a:  1  P1.A2 

6010  221 

216  42  •  1U0PI>A2 

221  42  •  424Pt 

42  •  42*TuOP1 
IP  1421  216.214,214 
216  42  •  424TU0PI 
214  PCIUAN 
CNO 
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INITIAL  DISTRIBUTION  LIST 


1.  Defense  Technical  Information  Center  2 

Cameron  Station 

Alexandria,  VA  22314 

2.  Library,  Code  0142  2 

Naval  Postgraduade  School 

Monterey,  CA  93943 

3.  Chairman  (Code  68Co)  S' 

Department  of  Oceanography 

Naval  Postgraduade  School 
Monterey,  CA  93943 

4.  Prof.  Stevens  P.  Tucker  ^ 

Depertment  of  Oceonography  (Code  68Tx) 

Naval  Postgraduade  School 
Monterey,  CA  93943 

5.  Director  Naval  Oceanography  Division  1 

Naval  Observatory 

34th  and  Massachusetts  Avenue  NW 
Washington,  DC  20390 

6 .  Commander  1 

Naval  Oceanography  Command 

NSTL  Station 

Bay  St.  Louis,  MS  39522 

7.  Commanding  Officer  1 

Naval  Oceanograph ic  Office 

NSTL  Station 

Bay  St.  Louis,  MS  39522 

8.  Commanding  Officer  1 

Fleet  Numerical  Oceanography  Center 

Monterey,  CA  93943 

9.  Commanding  Officer  1 

Naval  Ocean  Research  and  Development 

Activity 

NSTL  Station 

Bay  St.  Louis,  MS  39522 

10.  Commanding  Officer  ^ 

Naval  Environmental  Prediction 
Research  Facility 
Monterey,  CA  93943 
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11.  Chairman,  Oceanography  Department  1 

U.S.  Naval  Academy 

Annapolis,  MD  21402 

12.  Chief  of  Naval  Research  1 

800  N.  Quincy  Street 

Arlington,  VA  22217 

13.  Director  (Code  PPH)  I 

Defense  Mapping  Agency 

Bldg.  66,  U.S.  Naval  Observatory 
Washington,  DC  20305 

14.  Director  (Code  HO)  3 

Defense  Mapping  Agency  Hydrographic 
Topographic  Center 

6500  Brookes  Lane 
Washington,  DC  20315 

15.  Director  (Code  PSD-MC)  1 

Defense  Mapping  School 

Ft.  Belvoir,  VA  22060 

16.  Director,  Charting  and  Geodetic  1 

Services  (N/CG) 

National  Ocean  and  Atmospheric 

Administration 

Rockville,  MD  20852 

17.  Chief,  Program  Planning,  Liaison  1 

and  Training  (NC2) 

National  Oceanic  and  Atmospheric 

Administration 

Rockville,  MD  20852 

18.  Chief,  Nautical  Charting  Division  (N/CG2)  1 

National  Oceanic  and  Atmospheric 
Administration 

Rockville,  MD  20852 

19.  Chief,  Hydrographic  Surveys  Branch  1 

(N/CG24 ) 

National  Oceanic  and  Atmospheric 

Administration 

Rockville,  MD  20852 

20.  Director,  Pacific  Marine  Center  1 

(N/MOP) 

National  Ocean  Service,  NOAA 
1801  Fairview  Avenue  East 
Seattle,  WA  98102 
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21.  Director,  Atlantic  Marine  Center  1 

(N/MOA) 

National  Ocean  Service,  NOAA 
439  W.  York  Street 
Norfolk,  VA  23510 

22.  Cdr .  Kurt  J.  Schnebele,  NOAA  1 

Dept,  of  Oceanography  (Code  68SN) 

Naval  Postgraduate  School 
Monterey,  CA  93943 

23.  Capt.  Michael  lacona  (USNR  ret)  1 

2304  Roslyn  Ave . 

District  Hts . ,  MD  20747 

24.  Comandant  1 

G-NRN-1 

2100  SW  2nd  St.  SW 
Washington,  DC  20593-0001 

25.  Defense  Mapping  Agency  2 

Hydrographic/Topographic  Center 

Code  HCNN,  Attn:  L.  Phillips 
6500  Brookes  Lane 
Washington,  DC  20315 

26.  IHO/FIG  International  Advisory  Board  1 

International  Hydrographic  Bureau 

Avenue  President  J.F.  Kennedy 
Monte  Carlo,  Monaco 

27.  Hydrographer  of  the  Navy  1 

Ministry  of  Defense 

Hydrographic  Department 
Taunton,  Somerset 
TAI29N  England 

28.  Hydrographer,  Royal  Australian  Navy  1 

Hydrographic  Office 

Box  1332 

North  Sydney  2060 

New  South  Wales,  Australia 

29.  Dominion  Hydrographer  1 

Canadian  Hydrographic  Service 

615  Booth  Street 
Ottawa,  Ontario 
Canada  K1AOE6 

30.  Kapala  Jawatan  Hidro  Oseanografi  1 

Jalan  Gunung  Sa)iari  -  87 

Jakarta,  Indonesia 
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31.  Dlreccion  General  De  Oceanografia 
Departamento  De  Hidrograf icia 
Ave.  COYOACAN  131 
Mexico  12,  D.F. 

32.  Director  1 

Instituto  Hidrografico 

Rua  Das  Trinas,  49 
Lisboa,  Portugal 

33.  Civil  Engineering  Department  1 

Ports  and  Lighthouses  Administration 

Ras  Eltin  -  Gate  1 
Alexandria,  Egypt 

34.  Hydrographic  Service  1 

Athens  BST  902,  Greece 

35.  Chief,  Hydrographic  Surveys  Branch  1 

Atlantic  Marine  Center,  NOAA 

439  West  York  St. 

Norfolk,  VA  23510 

36.  Commanding  Officer  1 

Naval  Western  Oceanography  Center 

Box  113 

Pearl  Harbor,  HI  96860 

37.  Commanding  Officer  1 

Naval  Oceanography  Command  Center,  Rota 

Box  31 

FPO  San  Francisco,  CA  96540 

38.  Library  1 

Moss  Landing  Marine  Lab 

California  State  Colleges 

Sandholdt  Road 

Moss  Landing,  CA  95039 

39.  Library  2 

School  of  Oceanography 

Oregon  State  University 
Corvallis,  OR  97331 

40.  Library  1 

Dept,  of  Oceanography 

University  of  Washington 
Seattle,  WA  98105 

41.  ARNAV  Systems  Inc.  (attn:  Mr.  Walter  Dean)  1 

4740  Ridge  Dr .  NE 

Salem,  OR  97303-0012 
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42.  Regional  Director 

U.S.  Dept,  of  the  Interior 
Minerals  Management  Service 
1340  W.  6th  St. 

Los  Angeles,  CA  90017 
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